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COUNTER  M0DEU3  AMD  APPLICATIONS  TO  DETECTION  PROBLEMS 

by 

Hans  Zweig 


INTRODUCTION  AND  SUMMARY 

Work  in  several  aspects  of  detection  theory  has  been  in  progress 
over  the  last  two  decades.  Electronic  engineers  dealing  with  communi¬ 
cation  theory  have  looked  into  detection  problems  which  develop  from 
their  studies.  Another  direction  of  work  in  detection  theory  stems 
from  the  work  of  physicists  in  the  area  of  radiation  and  infrared  de¬ 
tection.  Detection  theory  is  a  term  also  applied  in  research  problems 
and  in  other  contexts.  The  concepts  and  methods  of  probability  and 
statistics  play  important  roles  in  these  studies  and  vary  with  the 
character  of  the  problems.  Communication  theory  is  most  closely  allied 
with  the  branch  of  statistics  dealing  with  hypothesis  testing  and  the 
work  of  transposing  this  statistical  methodology  into  the  field  of  com¬ 
munication  engineering  receives  much  attention.  The  second  study  men¬ 
tioned  above  is  concerned  more  with  estimation  problems  and  has  not  had 
such  a  straightforward  development  as  the  first.  It  has  not  been  clear 
which  physical  and  statistical  concepts  are  useful  nor  how  detectors 
(l.e.>  radiation  detectors)  should  be  evaluated.  Here  the  major  problem 
is  one  of  determining  how  actual  detector  devices  fall  short  of  being 
ideal  in  some  sense.  In  this  case  more  than  in  communication  theory 
there  is  a  need  to  simulate  mathematically  or  probabilistically  the 
working  of  the  detectors  in  order  to  understand  and  evaluate  the  detec¬ 
tion  process. 


Radiation  la  oftan  thought  of  aa  a  Rolaaon  stream  of  discrete  quanta 
or  events,  or  a  combination  of  Poisson  streams,  and  radiation  detectors 
can  therefore  be  thought  of  as  estimators  of  the  parameter  or  parameters 
of  these  streams.  One  criterion  for  Judging  the  performance  of  radiation 
detectors  which  appears  particularly  useful  Is  that  of  "quantum  effi¬ 
ciency"  .  This  expression  has  been  used  in  a  variety  of  ways.  Like  most 
Indices  it  has  grown  out  of  Intuitive  notions  and  thus  some  Initial  de¬ 
finitions  appropriate  for  the  original  motivation  must  be  used  carefully 
in  other  settings.  In  Chapter  I  we  give  a  specific  definition  which  is 
used  throughout  the  paper. 

The  history  of  this  expression  is  interesting.  Originally  it  was 
used  simply  to  designate,  for  a  given  device,  the  ratio  of  the  number 
of  output  events,  M,  to  the  number  of  input  events,  Nj  a  usage  which 
can  be  shown  to  be  quite  inappropriate  in  a  large  number  of  situations 
(see  Chapter  IV).  The  modern  version  of  this  concept  was  first  formu¬ 
lated  by  Rose  (ll)  and  later  reintroduced  under  the  designations  "detec¬ 
tive  quantum  efficiency"  by  Jones  (5,6)  and  "quantum  efficiency"  by 
Fellgett , ( k) .  These  writers  defined  the  notion  as  a  ratio  of  signal-to- 
nolse-ratios  of  the  "ideal  detector"  to  that  of  the  actual  detector. 
Another  definition  was  proposed  by  Zveig,  Higgins,  and  MacAdam  (20). 

It  is  the  ratio,  lying  between  zero  and  one,  of  variances  of  estimates 
of  intensity  (e.g.  mean  of  a  Bcisson  process)  for  an  ideal  detector 
(numerator)  and  for  an  actual  device  (denominator).  They  also  demon¬ 
strated  that  this  definition  is  equivalent  to  that  used  by  the  earlier 
writers  when  it  is  applied  to  their  analyses.  Using  this  concept  of 
quantum  efficiency,  Zveig  (18,  19)  examined  a  simplified  model  of  the 
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Ale  emcaso  to.  atadrbagfP«saiaa4lialtatloas  of  tfirxoasaefcle'  • 
detector*. 

The  purpose  of  the  present  work  Is  to  derive  the  quantum  efficiency 
of  a  variety  of  detectors  for  which  mathematical  models  already  exist 
and  to  develop  nev  models  for  some  physical  situations  which  can  be 
examined  in  this  light. 

The  plan  of  this  work  is  as  follows:  In  Chapter  I,  (a)  a  summary 
of  the  models  which  will  be  analyzed  is  given,  (b)  the  applications  of 
these  models  to  specific  detection  problems  is  Indicated,  and  (c)  the 
methods  of  estimating  the  parameter  of  Interest  and  the  efficiency  of 
the  estimates  is  defined.  Chapter  II  contains  a  summary  of  known  dead 
time  counter  models  and  evaluates  their  quantum  efficiency.  In  Chapter 
III  a  nev  class  of  counters,  called  coincidence  counters,  is  Introduced 
and  evaluated.  Chapter  IV  contains  a  discussion  of  photographic  detec¬ 
tors.  Hie  dead  time  and  coincidence  counter  models  are  used  here  to 
simulate  the  photographic  effect  known  technically  as  "reciprocity 
failure",  and  the  quantum  efficiency  of  photographic  detectors  both  with 
and  without  reciprocity  failure  is  obtained.  Chapter  V  contains  two 
other  applications  of  the  counter  models.  In  this  chapter  it  is  shown 
how  combinations  of  counters  can  be  used  to  solve  more  complex  deteotlon 
problems.  One  application  involves  using  dead  time  and  coincidence 
counters  in  aeries.  Ibis  combination  is  useful  in  solving  a  radar 
ranging  problem  when  the  detection  problem  consists  of  determining  the 
presence  of  a  pulse  pair  in  a  Poisson  nolee  stream.  Another  application 
involves  using  these  counters  in  parallel*  1°  this  case  the  application 
of  interest  la  that  of  determining  the  component  parameters  of  a  compound 
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1.  •  Model* 

Various  types  of  counter  models  have  been  Introduced  and  discussed 
In  the  literature  of  the  last  tventy  years .  Type  I  and  Type  II  counters 
have  received  the  most  extensive  treatment.  Type  I  counters  are  charac¬ 
terized  by  the  fact  that  a  registered  event  produces  a  dead  time  during 
which  no  further  events  can  be  registered.  Type  II  counters  are  charac¬ 
terized  by  the  fact  that  incoming  events  produce  dead  time  during  which 
further  incoming  events  cannot  be  registered,  although  these  events 
are  capable  of  prolonging  the  dead  time,  that  is,  paralyzing  the  counter. 
Precise  definitions  of  the  counter  models  will  be  found  in  the  sections 
in  which  these  models  are  discussed  and  evaluated.  Modifications  of 
Type  I  and  Type  U  counters  have  also  been  introduced  and  studied  in  the 
literature . 

Although  the  problem  of  dealing  with  transformations  of  Poisson 
processes  goes  back  at  least  as  far  as  Bortkiewiez ' s  paper  of  1913  (l), 
the  earliest  formulation  of  what  has  come  to  be  called  the  Type  II 
counter  seems  to  be  that  of  Levert  and  Scheen  in  19^3  (7) •  The  distinc¬ 
tion  between  Type  I  and  Type  II  counters  was  made  somewhat  earlier  (1937) 
by  Ruark  and  Brammer  (12).  Malmquist  (8)  considered  the  case  of  a  Type 
I  counter  operating  on  a  process  with  arbitrary  inter-arrival  time  dis¬ 
tribution  f(t),  thus  generalizing  the  Poisson  input  inter-arrival  that 
distribution  F(t)  «  l-e“XT  rather  than  the  counter  mechanism.  Feller 
(3).  confining  himself  to  the  classical  Type  I  and  Type  II  models,  has 


jr*»*»ted  tin  «m  of  straightforward  probabilistic  arguments  and  Laplace 
irtasforms  to  determine  the  I&place  transform  of  the  inter-arrival  time 
distribution  of  registered  events.  Takacs  (l 6),  using  a  similar  approach, 
generalised  the  counter  model  so  that  each  incoming  event  produces  a 
dead  time  Interval  t  vith  probability  p  and  another  event  occurring 
during  this  dead  time,  if  one  is  created,  will  not  be  registered  although 
it  may  also  produce  another  dead  time  with  the  same  probability.  The 
case  p  ■  0  produces  the  Type  1  counter,  while  p  «  1  produces  the 
Type  II  model.  In  this  paper,  we  shall  evaluate  the  quantum  efficiency 
of  Type  I  and  Type  II  counters.  These  provide  Interesting  comparisons 
with  the  coincidence  counter  models  which  will  be  developed,  as  well  as 
with  the  photographic  type  of  detectors. 

Takacs1  model  is  interesting  in  that  it  provides  a  bridge  between 
the  Type  I  and  Type  II  counters.  However,  in  doing  this  it  also  blurs 
the  distinctive  features  of  these  counters  and  it  is  these  that  we  want 
to  exploit.  Takacs  has  also  dealt  with  a  coincidence  problem  which  is 
concerned  with  the  probability  of  finding  several  dead  time  counters 
simultaneously  locked.  This  is  a  different  type  of  coincidence  situa¬ 
tion  than  the  one  we  now  present. 

The  simple  coincidence  counter  model  we  now  introduce  is  Intended 
to  work  as  follows:  An  incoming  event,  say  a  photon,  is  not  directly 
recorded  by  the  detector.  Instead  it  changes  the  state  of  the  detector 
say  from  SQ  to  S^.  The  stats  8^  is  maintained  for  a  fixed  time 
following  the  incidence  of  the  photon.  If  during  this  time  interval  no 
further  arrival  takes  place,  the  detector  reverts  to  state  SQ  and  rjo 

:  --'r,  ,  ;  t  .  '  ‘  ,  X  :  '  j'  v-  i  ^ 

output  count  is  recorded.  If,  however,  while  the  detector  is  in  state 
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t  *  yt  as  output  avant  would  be  registered  at  tines  t  +  2. It  and  at 
t  ♦  5t|  or  in  thii  casa  five  input  events  have  produced  two  output  event*. 

Xhe  R-fold  coincidence  counter  la  a  generalisation  of  the  simple 
coincidence  counter.  Here,  incoming  events  cause  changes  in  the  state 
of  the  detector  from  to  S1+1  provided  the  detector  is  in  state 
Sj  (i  »  0,1, ...,R-l)  at  the  time  of  the  arrival  of  the  input  event. 

For  i  »  1,2, ...,R-1  the  counter  remains  in  state  i  for  a  time  t. 
which  is  now  taken  to  he  a  random  variable  with  distribution  indepen¬ 
dent  of  i .  If  no  further  event  occurs  during  time  following  the 
last  input  the  counter  reverts  to  state  SQ  without  registering  a  count. 
Only  if  state  Sg  is  reached  is  a  count  registered  with  the  counter 
imaediately  reverting  to  SQ. 

The  photographic  detector  is  a  counting  device  of  an  entirely  dif¬ 
ferent  type  but  the  counter  models  nevertheless  play  a  part  in  the  ana¬ 
lysis  of  this  type  of  detector.  Basically  the  photographic  detector 
(the  emulsion  or  photographic  plate)  consists  of  an  ensemble  of  many 
counters  (the  individual  photographic  grains),  each  of  which  is  able 
to  register  a  count  only  once.  The  fractional  ntaaber  of  grains  that 
develop  is  then  a  measure  directly  related  to  the  average  mother  of 
Photon*,  incident  on  eenh  grain  during  the  time  of  eaqposure.  The  slat- . 
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and  speed,  i*e.,  each  grain  requires  B  photons  to  make  it  developable 

■  '  .  '  !  ?  1  '  *  *'  1  JV,  y  ■ 

(that  is,  to  register  a  count).  Aside  from  vmristiene  ln  graia  eise  ' 
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which  will  not  be  considered,  actual  photographic  detectors  consist  of 
an  ensemble  of  grains  of  varying  speed  so  that  R,  the  number  of  photons 
required  for  registering  an  output  event,  is  a  random  variable  with 
distribution  dependent  on  the  particular  photographic  material. 

2.  Applications 

There  is,  a  photographic  effect  which  can  be  explained  in  terms  of 
the  dead  time  and  coincidence  counter  models  introduced  above.  This 
effect  is  known  technically  as  reciprocity  failure  and  manifests  itself 
in  two  different  ways,  (a)  At  high  intensities  X.  and  short  exposure 
times  T  the  photographic  detector  registers  fewer  events  for  a  given 
average  number  of  photons  per  grain  XT  than  when  X  is  somewhat 
smaller  (but  XT  is  fixed).  This  can  be  explained  by  postulating  a 
dead  time  t  following  each  incident  photon  so  that  only  those  photons 
contribute  to  produce  the  necessary  number  R  of  input  events  required 
for  developability  of  the  grain  which  arrive  at  least  t  units  apart, 
(b)  Low  intensity  reciprocity  failure  denotes  the  condition  in  which 
the  response,  for  a  fixed  average  input  XT,  is  again  lowered,  this 
time  when  the  intensity  of  radiation  X  is  low;  that  is,  the  photons 
arriving  at  a  grain  tend  to  be  spaced  far  apart.  By  assuming  the  type 
of  mechanism  discussed  in  the  case  of  the  coincidence  counter,  this 
type  of  reciprocity  failure  can  also  be  explained.  The  dead  time 
counter  models  and  coincidence  time  counter  models  postulated  for  photo¬ 
graphic  detectors  can  be  identified  with  the  times  required  for  the 
photochemical  creation  of  electron  hole  pairs  and  the  recombination  for 
such  pairs  respectively. 
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Apart  from  thie  application  of  the  counter  models,  other  possibi¬ 
lities  suggest  themselves.  The  use  of  a  Type  I  and  a  simple  coincidence 
counter  in  series  enables  one  to  construct  a  device  for  detecting  the 
presence  of  a  pulse  pair  of  known  spacing  embedded  in  a  Poisson  noise 
stream,  nils  combination  of  counters  suggests  itself  from  range  deter¬ 
mination  problems  in  space  systems.  Another  possibility  would  be  to 
use  these  counters  in  parallel .  If  the  Poisson  stream  consists  of  a 
variety  of  particles  of  different  energy  E,  and  the  coincidence  counter 
can  be  constructed  so  that  particles  of  energy  E  <  Eq  produce  only 
one  transition  (SQ  -*  S^)  whereas  particles  with  energy  E  >  E^ 
produce  two  or  more,  then  the  coincidence  counter  will  respond  primarily 
to  the  high  energy  particles.  A  dead  time  counter,  on  the  other  hand, 
will  not  differentiate  between  particles  of  different  energy.  By  varying 
the  level  Eq  it  is  thus  possible  to  estimate  the  energy  distribution 
of  the  Poisson  stream.  Since  energy  is  directly  proportional  to  fre¬ 
quency  what  we  have  thus  described  is  the  discrete  counterpart  of  a 
spectrum  analyzer.  This  analysis  may  be  useful  in  the  study  of  separa¬ 
tion  of  energy  levels  for  high  energy  particles. 

3.  Methods 

The  method  for  obtaining  estimates  of  the  Poisson  stream  parameter 
is  the  following.  We  take  as  our  datum  the  quantity  M(t)  which  is 
the  observed  total  number  of  output  counts  from  the  detector  in  time  T. 

From  M(T)  we  need  an  estimate,  say  f[M(t)],  of  the  expected 
number  of  counts,  given  X,  in  time  T,  E{M^(T)).  Asymptotically, 
i.e.,  for  large  T,  we  know  by  the  law  of  large  numbers  that  M(T)  will 
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approach  E(M^(T) )  provided  the  variance  of  M(T)  is  finite.  Now 
since  M(T)  Is  an  asymptotically  unbiased  estimate  of  E{M^(T))  we 
choose  f [M(T)]  ■  M(T).  Thus  M(T)  is  our  estimate  of  the  expected 
number  of  output  counts.  The  mathematical  model  will  provide  a  relation 
between  the  parameter  X  of  the  Poisson  stream  and  the  expected  number 
of  output  counts  in  a  fixed  time  of  observation  T,  namely 

A 

r(X)  =  E(M^Ct)).  We  then  obtain  an  estimate  of  X,  call  it  X,  by 
equating  M(T)  =  r(X)  and  solving  for  X.  This  yields  a  unique, 
asymptotically  unbiased  estimate  provided  the  expression  r(X)  is 
continuous  and  monotonic  in  X.  In  the  case  of  counters  the  expression 
for  E(M^(T) }  is  derived  from  the  Laplace  transform  of  the  inter¬ 
arrival  time  random  variable  using  results  from  renewal  theory.  In  the 
case  of  the  photographic  detector  the  expression  for  E{M^(T) )  1b  de¬ 
rived  by  direct  argument. 

The  ideal  detector  is  one  whose  estimate  of  X  call  it  X*  has 
the  lowest  possible  variance  in  the  class  of  unbiased  detectors.  When 
we  are  dealing  with  Poisson  inputs  we  can  achieve  a  variance  whose 
lower  bound  is  given  by  the  Cramer-Rao  Theorem.  In  this  case  we  know 
that  a  uniform  minimum  variance  estimator  of  XT  ,  when  the  process  is 
observed  for  a  total  time  T,  is  N(T),  the  total  number  of  input 
events.  This  estimate  has  variance  XT  and  hence  the  variance  of  the 
best  estimate  of  X,  i.e.,  X*  ,  is  obtained  from 

(1.3.1)  var{X*}  =  varlnjr^}  =  ^  var{N(T)}  =  7-  . 

Having  settled  on  a  method  for  estimating  X  and  an  upper  bound 
on  the  precision  (lower  bound  on  variance)  attainable  from  any  detector 
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we  can  define  a  measure  of  efficiency  vhich  will  allow  us  to  compare 
different  detectors . 

We  define  the  quantum  efficiency  of  a  detector  as  the  ratio  of  the 

A 

variance  of  the  estimate  X*  to  the  variance  of  the  estimate  X  pro¬ 
duced  by  the  given  detector.  This  efficiency  is  a  function  of  both  X 
and  the  length  of  the  time  interval  of  observation  T,  (and  other  para¬ 
meters  depending  on  the  model)  although  in  many  cases  it  will  be  inde¬ 
pendent  of  T . 

We  shall  be  concerned  here  with  asymptotic  results  and  we  Bhall 
obtain  the  asymptotic  quantum  efficiency  in  the  following  way:  We 
obtain  asymptotic  results  for  E{.'M^(T))  -  r(X)  and  then  obtain  the 

A 

asymptotic  estimate  X  by  solving  the  equation  M(T)  =  r(X).  To  obtain 
the  asymptotic  variance  of  this  estimate  we  expand  M(T)  =  r(X),  in  a 
Taylor  Series  about  X,  the  true  value  of  the  parameter, 

(1.3.2)  M(T)  =  r(X)  =  r(X)  +  (X-X)  +  . . 


or 


(1.3.3) 


M(t)  -  r(X)  = 


d  E{M^(T)} 


dX 


(X-X) 


Squaring  both  sides  of  (1.3-3)  and  taking  expectations  (ignoring 
higher-order  terms)  we  find  the  variance  of  M(T)  is  related  to  the 
variance  of  the  estimate  of  X  by  the  equation 

(1.3. 4)  var(M(T))  =  var{X) 
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or 


(1.3.5) 


*  d  EfM^T))  -2 


var{X)  =  [• 


dX 


]  •  var{M(T))  . 


Since  the  ideal  detector  estimates  X  by  X*  which  we  know  has 
variance  X/t  we  note  that  the  asymptotic  efficiency,  is  given  by 


(1.3.6) 


c  _  var(X»  1  N(T) ) 
*  var{£  j  M(T) ) 


or 


(1.3.7) 


d  E{M,(T)}  2 

r - 1 - j 

c  =  b  .  dX  * 

*  var{M(T)) 


The  discussion  in  this  study  is  concerned  largely  with  such  de¬ 
vices  as  the  radiation  detectors  and  counting  tubes  used  and  studied 
by  physicists,  and  for  these  the  measure  £  is  extremely  useful. 

Other  measures  of  performance  would  need  to  be  devised  for  detectors 
whose  function  it  is  to  decide  between  alternative  states  of  nature. 
For  example,  in  the  case  of  two  such  states,  a  criterion  of  efficiency 
might  involve  not  only  the  probability  of  a  correct  decision,  but  the 
time  involved  in  reaching  a  decision.  This  would  be  especially  true 
in  monitoring  fallout,  where  the  object  is  to  decide  whether  the  level 
is  or  is  not  critical.  A  simpler  criterion  for  evaluating  a  detector 
which  must  decide  which  of  two  possible  states  prevails  would  be  the 
total  probability  of  error  over  a  given  time  span.  We  shall  evaluate 
a  particularly  interesting  detector  with  respect  to  this  criterion  in 
Chapter  V. 


12 


CHAPTER  II 


THE  DEAD  TIME  COUNTERS 

1.  Results  From  Renewal  Theory 

We  shall  follow  Feller  (3)  in  deriving  the  Laplace  transform  of 
the  inter-arrival  time  random  variable  for  the  Type  I  and  Type  II 
counters,  since  the  type  of  argument  used,  particularly  for  the  Type  II 
counters,  lends  itself  to  deriving  a  corresponding  result  for  a  simple 
coincidence  counter  mechanism.  However  before  deriving  these  results, 
we  will  look  into  some  aspects  of  renewal  theory  which  will  be  helpful 
in  arriving  at  the  efficiency  of  the  detectors  under  consideration. 

What  we  require  is  a  relation  between  the  mean  and  the  variance  of 
the  inter-arrival  distribution  and  the  mean  and  variance  of  the  number 
of  events  that  will  occur  in  a  given  time.  These  results  are  discussed 
by  Smith  (lk),  and  an  elementary  treatment  is  given  by  Saaty  (13).  Let 
Xi  be  the  random  variable  which  denotes  the  elapsed  time  between  suc¬ 
cessive  events  (renewals),  all  of  which,  except  perhaps  the  first,  are 
assumed  identically  distributed.  Then,  we  can  write 

(2.1.1)  T„-Xl+"-*X» 

to  denote  the  waiting  time  random  variable  up  to,  and  including,  the 
n  event.  This  waiting  time  random  variable  has  a  distribution  re¬ 
lated  to  the  distribution  function  F(t)  of  X^,  which  is  given  by 

(2.1.2)  F^(t)  =  f  F^n_1^(t-x)  cLF(x) 

Jo 

(n) 

where  F'  ;(t)  denotes  the  n-fold  convolution  of  F(t)  with  itself. 
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The  are  also  assumed  to  be  Independently  distributed. 

Let  exactly  M  events  take  place  In  a  given  time  T.  Label  this 
M(T).  Tines  probability  of  this  happening  is 

(2.1.3)  PM  =  F^(T)  -  j'(Mfl)(T) 

since  this  is  the  probability  that  the  event  but  not  the  (WH)8^ 
event  takes  place  prior  to  T.  From  we  can  then  derive  the  mean 
and  variance  of  the  number  of  events  in  (0,T).  Letting  E{M(T))  and 
var{M(T)}  denote  these  parameters,  and  |i  and  a  the  mean  and 
variance  of  the  distribution  F(T),  we  have,  from  the  key  renewal 
theorem  of  Smith  (19)  that,  as  T  -» ao 

(2.1. 4)  E(M(T))»^ 
and 

2 

(2.1.5)  var{M(T)}  =-2—  . 

H 

These  results  can  be  viewed  in  the  following  non-rigorous  manner. 
Heurlstically  we  may  write 

(2.1.6)  E{X1  +  ...  +  XmW)  =  n  E{M(T) }  , 

that  is,  the  expected  time  for  M  events  should  approximate  the  expected 
time  for  one  event  multiplied  by  the  expected  number  of  events.  Likewise 
we  may  also  write 

var{X1  +  ...  +  X^}  =  o2E(M(T)} 

remembering  this  is  an  approximation  based  on  the  same  argument  since 
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the  variance  of  a  sum  of  independent  random  variables  Is  the  sum  of  the 
variables.  Thus  ve  may  also  vrlte  In  the  same  approximate  vein 

(2.1.7)  £[3^  +  ...  +  Xn(T)  -  u  M(T)]2  -  a2  B(M(T)) 

Furthermore,  for  large  T  we  expect 

(2.1.8)  ^  +  ...  +  XM(T)  -  T  . 

The  actual  value  is  of  the  order  T  -  \i,  (recall  n  is  the  mean  waiting 
time  for  an  event  and  thus  for  the  first  event)  which  for  large  M  can 
be  approximated  by  T,  hence  using  (2.1.6) 

(2.1.9)  E{M(T))  -»2 

and  substituting  in  (2.1.7)  we  get 

'  2 

(2.1.10)  E{T  -  u  M(T)}2  -  ip  . 

Multiply  both  sides  by  1/^2  obtaining 

(2.1.11)  E[M(T)  -  J}2  «  £_| 

and  since  ^  can  be  replaced  by  E{M(T)J,  we  get 

ct^T 

(2.1.12)  var{M(T)}  -  -2-i  . 

With  these  results  in  mind,  we  return  to  the  evaluation  of  counters. 

2.  The  Type  I  Counter 

The  operation  of  a  Type  I  counter  can  be  described  briefly  as 
follows: 
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A  counter  which  immediately  registers  the 
first  input  event  and  then  only  those  input 
events,  which  are  not  preceded  during  a 
fixed  time  interval  hy  a  registered  (output ) 
event . 


The  results  for  the  Type  1  counter  are  those  which  are  most  easily 

obtained.  The  dead  time  is  assumed  to  he  a  constant  t  and  the  counter 

operates  on  a  Poisson  process  which  has  for  inter-arrival  distribution 

-X.t 

of  input  events  the  negative  exponential  distribution  1-e  .  Thus 

the  inter-arrival  time  of  registered  events  have  a  distribution  (after 
occurrence  of  the  first  event) 


(2.2.1) 


ro(t) 


-X(t-x) 


t  <  T 
t  >  X 


FQ(t)  is  used  throughout  to  describe  the  output  inter-arrival  time 
distribution,  F(t)  has  been  and  will  continue  to  be  used  for  the 
input  inter-arrival  time  distribution. 

The  Laplace  transform  of  t  is 


(2.2.2) 


Since 

(2.2.3) 

and 


=  J  t  dFQ(t) 


d  cpo(s)  ( 

8=0 
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(2.2.4) 


2.2 
o+4 
o  o 


r"  2  d2  *>> 

[  t2  dP  (t) - ^5 — 

Jo 


dt 


8»0 


we  obtain  the  mean  and  variance  of  the  inter-arrival  distribution  of 
counted  events: 


%  = 


Vr+i 

x 


(2.2.5) 

and 

(2.2.6)  a2  =  -i 

0  X* 

Employing  (2.1.4)  and  (2.1.12),  we  get 

(2.2.7)  E(M(T) )  =-5~ 


and 

(2.2.8)  var(M(T))  =  -  H  -  . 

(Xt+1)5 

Before  we  compute  the  quantum  efficiency,  it  is  interesting  to 
note  that  for  large  values  of  X,  that  is,  strong  incident  radiation, 
we  get 


(2.2.9) 


E{M(T) )  -»  T/t 


and 

(2.2.10)  var(M(T)J  -+ 0  . 

This  result  can  be  viewed  in  the  following  way.  In  the  limit  as  X  -» co 
we  get  counts  periodically  and  immediately  following  the  end  of  the  last 
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dead  time,  so  that  the  number  of  counts  becomes  strictly  proportional 
to  the  observation  time,  independent  of  X,  and  thus  the  variance  is 


zero. 


The  quantum  efficiency  is  obtained  from  the  expression  (1.3*7) 


(2.2.11) 


C  X  rdgfM(T))  I2 _ 1 

£  "  7  L  dX  J  varfM 


var 


nm 


which,  in  this  case,  turns  out  to  be 

„2 


(2.2.12) 


c  _  X  T 

yp  ~  rjj 


(Xt+1) 


(xrjij; 

5  XT 


or 


(2.2.13) 


£  =  Xt+1  * 


We  see  that  for  non- zero  t  the  quantum  efficiency  tends  monotonically 
to  zero  as  X  ->oo;  that  is,  for  strong  radiation,  the  Type  I  counter 
is  not  very  efficient  as  an  estimator  of  X.  The  number  of  counts 
increases  with  the  time  of  observation,  but  tends,  for  large  X,  to 
become  independent  of  the  radiation  strength  and  only  a  function  of  the 
dead  time.  Hie  efficiency  depends  on  the  arrival  rate  X  and  is  in¬ 
dependent  of  the  time  of  observation  T  or  the  total  number  of  incident 
events . 


3.  The  Type  II  Counter 


Hie  operation  on  the  Type  II  counter  can  be  stated  as  follows: 

A  counter  which  immediately  registers  those, 
and  only  those  input  events,  which  are  not 
preceded  during  a  fixed  time  interval  by  an 
input  event. 
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We  follow  Feller's  (3)  argument  is  establishing  the  inter-arrival 
time  distribution  of  output  events  for  the  Type  II  counter.  It  should 
be  noted  Just  from  the  physical  situation  that  the  number  of  registered 
events  in  the  presence  of  strong  incident  radiation  tends  to  zero  for 
any  given  time  interval  T. 

Let 

(2.3.1)  P  =  e'XT 
and 

(2.3.2)  q  =  l-e-XT 

denote  respectively  the  probability  of  a  non-arrival  and  an  arrival  of 

an  incident  event  during  a  time  span  t.  Then  the  probability,  once  a 

counter  is  paralyzed,  that  V  events  will  prolong  the  dead  time  is 
v 

q  p.  Let  X^  denote  the  elapsed  time  between  input  event  numbers 
(i-l)  and  i.  Then  we  can  write 

(2.3.3)  T(v)  =  X±  +  •••  +  Xv  +  t 

for  the  total  time  that  the  counter  is  paralyzed  by  exactly  V  events. 

This  is  a  random  variable,  conditioned  on  v,  which  is  the  sum  of  V 

independent,  identically  distributed  random  variables.  It  follows 

••ST 

that  the  Laplace  transform  of  T(v)  is  the  product  of  e  and  the 
V  power  of  the  Laplace  transform  of  X^ .  Once  we  have  the  transform 
of  T(v),  we  can  multiply  qVp  and  sum  over  V  to  obtain  the  Laplace 
transform  of  the  unconditional  paralyzed  time  random  variable  which  we 
need  to  evaluate  the  Type  II  detector. 
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Hie  probability  that  an  event  vlll  occur  in  a  time  (0,t)  with 
t  <  t  given  that  it  occurs  in  a  time  interval  (0,t)  is 


(2.3.4) 


l-e-Xt  l-e-Xt 


and  this  is  the  conditional  probability  <  t),  for  i«l,  . . . ,  v 


and  gives  rise  to  the  Laplace  transform 


(2.3.5)  /’  e'"  7^7  "  -  i  •  £ 

Jo  1-e  * 

so  that  the  Laplace  transform  of  T(v)  is  simply 


(2.3.6) 


<1 


The  unconditional  Laplace  transform  is  obtained  by  summing 


(2.3.7) 


I  e- 
v=o 


Mi 

n  ' 


X_  ^^-(X+sJtjJ 


which  yields 


(2.3.8) 


U+s)e-(s+X,)T 
,  ,-  ( b+X)t 


s  +  X.e 


This  is  the  transform  of  the  paralyzed  time  distribution.  To 
obtain  the  transform  of  the  inter-arrival  time  of  registered  events,  we 
must  take  into  account  the  additional  time  from  the  end  of  paralysis 
to  the  arrival  of  another  incident  event.  The  latter  has  Laplace 

transform  t-~  so  that  we  obtain  finally  the  Laplace  transform  of  the 

A.+8 

inter-arrival  time  distribution  of  registered  events  as 
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(2.3.9) 


*o(b) 


JL_  .  *~(B+X)X 

X+B  b+X  e  (s+^'r 


■(b+X)t 


■  (s+X)t 


s  +  X  e 


The  next  step  is,  as  before,  to  obtain  the  mean  and  variance  of  the 
inter-arrival  time  distribution.  It  turns  out,  both  in  this  case  as  well 
as  others,  that  it  is  much  less  tedious  to  obtain  the  derivatives  of 
rather  than  tp(s) .  Noting  that 


(2.3.10) 


and  that  <p(0)  =  1  we  have 
(2.3.11) 

Also 


1*  >»'(0)  Ijjfej)  • 


(2.5.12)  -  -  <P"(s)  tsTTT)2  - 


ds 


l9TsTJ 


(<p(b)]' 


so  that  since 


(2.3.13) 


we  have 


+  a2  .»  [q>(0)] 

ds 


°2  -  -  'ifer'} 


(2.3.14) 

For  the  Type  II  counter  we  have 

(2.3.15) 


1  s  _ (X+b)t., 

TOT  ‘  *  e  41 
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so  that 


(2.3.16) 

and 


1  (X+s)t  st  (X.+s)t 
X  8  +  X  e 


(2.3.17) 


_af  .  1  ■  2t  (x+b)t  erf  (x+s)t 

d.2  x  1 


Consequently 


(2.3.18) 
and 

(2.3.19) 


Hence 


d_ 

ds 


1 

X 


XT 


e 


d2  r  1 
ds2 


] 


2t  Xt 

T  e 


(2.3.20) 


and 


(2.3.21) 


Again  using  (2.1.4)  and  (2.1.12)  we  obtain  the  mean  and  variance  of 
the  number  of  counts  in  an  interval  [0,T].  These  are 


(2.3.22) 


E{M(T) }  -  XT  e"XT 


and 

(2.3.23)  var (M(T) }  -  XT  e"2XT  [eXT  -  2Tt]  . 
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From  these  quantities  we  proceed  to  deduce  the  efficiency  of  the 
Type  II  counter.  Using  (2.2.11),  we  obtain  readily 

(2.3.24)  £  -  £ 

*  (eXT-2XT) 

The  character  of  this  result  suggests  that  we  look  at  the  Type  II 
counter  a  little  more  carefully.  Formula  (2.3.24)  asserts  that  the 
efficiency  of  this  type  of  counter  is  unity  when  X  is  zero,  decreases 
monotonically  to  zero  at  X  =  ^  and  then  becomes  positive  again  although 
it  tends  to  zero  for  X  =  oo  . 

This  type  of  behavior  is  not  quite  so  surprising  if  we  take  another 
look  at  the  distribution  of  the  output  counts,  M(T).  From  equation 
(2.3.22)  we  see  that  the  response  function  r(x)  *  E(M^(T)}  is  not 
monotonic  in  X.  As  the  radiation  strength  increases,  the  expected 
number  of  counts  increases  only  until  X  =  ^  .  For  X  >  ^  the  expected 
number  of  counts  actually  decreases  with  increasing  X,  since  the  amount 
of  time  that  the  counter  is  paralyzed  now  becomes  very  large.  A  con¬ 
sequence  of  this  is  that  the  function 

(2.3.25)  r(X)  =  XT  e"XT 

has  no  unique  inverse.  The  only  way  to  infer  a  unique  value  of  X 
from  a  single  observation  M(T)  is  to  have  additional  information 
about  X  which  guarantees  that  either  X  is  less  them  —  or  X  is 
larger  than  ~  .  In  the  absence  of  such  information  it  is  necessary  to 
have  an  auxiliary  experiment .  We  could,  for  example  introduce  ern 


absorption  device  into  the  Poisson  stream  which  reduces  its  intensity 
by  a  known  factor  Q.  The  reduced  stream  of  intensity  QX  will,  if 
X.  <  ^  )  produce  a  reduced  expected  number  of  counts  in  time  T.  If 
X  >  -  and  Q  is  near  unity  we  would  tend  to  get  a  larger  number  of 
output  events  from  the  reduced  stream.  This  auxiliary  experiment  would 
thus  allow  us  to  estimate,  uniquely,  the  parameter  X  of  the  Poisson 
stream.  We  see  from  this  that  it  should  not  usually  prove  to  difficult 
to  obtain  a  unique  estimate  of  X  from  two  observations  of  M(T) . 
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CHAPTER  III 


COINCIDENCE  COUNTERS 


1.  Introduction 

A  number  of  different  problems  in  a  variety  of  contexts  can  be  con¬ 
structed,  all  of  which  have  at  their  core  the  existence  of  an  integer 

valued  random  process,  X.(t),  which,  starting  from  the  state  S 

J  o 

strives  to  reach  some  preassigned  state  SQ.  The  process  can  increase 

K 

by  integer  amounts,  but  if  no  Jump  occurs  for  a  sufficiently  long  period 
of  time,  the  process  will  revert  to  its  initial  or  ground  state,  and 
must  begin  its  climb  towards  state  S_  all  over  again. 

We  can  imagine  a  psychologist  constructing  a  complex  learning 
or  maze  experiment  in  which  "success"  depends  on  the  performance  of 
R  successive  tasks,  and  if  any  task  is  not  performed  with  sufficient 
alacrity,  the  subject  must  begin  from  scratch.  Such  a  process  also 
occurs  in  the  study  of  photographic  detectors,  in  the  study  of  Lasers, 
and,  in  fact  in  a  great  variety  of  scientific  and  engineering  situations. 

The  fact  that  neither  Type  I  nor  Type  II  counters  (and  presumably 
intermediate  types)  are  very  good  detectors  for  large  values  of  X 
would  lead  one  to  suspect  that  in  those  cases  where  an  extreme  range 
of  parameter  values  X  is  a  priori  likely  to  occur,  different  counters 
ought  to  be  employed.  As  the  discussion  of  coincidence  counters  will 
show,  this  type  of  device  should  prove  of  some  value  in  such  cases  and 
moreover  it  provides  a  new  model  for  various  physical  situations  such 
as  those  mentioned  above. 
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2.  Ihe  Simple  Coincidence  Counter 


As  model  for  a  simple  coincidence  counter,  we  shall  consider  the 
following  mechanism:  An  Incoming  event  Impinges  upon  the  counter  which 
is  so  constructed  as  to  be  able  to  hold  this  event  In  memory  for  a 
time  t  (a  fixed  constant).  If  another  event  occurs  during  the  time 
that  the  first  event  is  in  memory,  the  counter  registers  an  output 
event  and  resets  its  memory  to  zero  content.  If  no  second  event  occurs 
during  time  t,  the  first  event  is  lost,  that  is,  produces  no  count. 

The  formal  definition  of  a  simple  coincidence  counter  can  be 
stated  as: 

A  counter  which  immediately  registers  an 
output  count  whenever  an  input  event  is 
preceded  during  a  fixed  time  interval  t 
by  a  nonregistered  input  event. 

Such  a  situation  may  be  obtained,  for  example,  on  a  photographic 
emulsion  in  which  the  first  photon  incident  on  a  silver  halide  grain 
creates  an  electron-hold  pair  which  has  a  lifetime  t  before  recom¬ 
bination  takes  place.  If,  however,  another  pair  is  formed  shortly 
after  the  first,  a  latent  image  speck,  that  is,  a  developable  silver 
nucleus,  or,  at  any  rate,  a  stable,  sublatent  image  speck  (i.e.,  a 
stable  speck  but  one  too  small  to  be  developed)  may  be  formed.  Another 
conceivable  situation  in  which  a  coincidence  counter  model  might  be 
useful  is  that  of  light  amplification  by  stimulated  emission  of  radi¬ 
ation  -  the  so-called  LASER.  In  this  device,  photons  are  required  to 
be  stored  until  a  sufficient  number  are  accumulated  which  are  then  re¬ 
leased  simultaneously  to  produce  an  intense  monochromatic,  coherent 
light  pulse.  If  an  insufficient  number  are  captured  during  the  holding 
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time  of  a  ruby  crystal  or  other  lasing  material  the  photons  which  arrive 
during  this  time  do  not  produce  a  light  pulse  and  are  wasted. 

It  should  be  noted  that  the  coincidence  model  suggested  above  dif¬ 
fers  from  the  general  counter  model  proposed  by  Takacs  in  that  the  oc¬ 
currence  of  an  event  within  t  units  of  an  earlier  one  produces  a  re¬ 
generation  point  (in  effect,  a  new  origin)  for  the  output  stochastic 
process  in  our  case,  but  not  in  Takacs'  case.  In  his  process,  an  output 
event  could  occur  after  the  arrival  of  both  the  second  and  the  third 
input  event.  In  our  model,  that  would  be  considered  physically  im¬ 
possible.  To  get  two  output  counts,  at  least  four  input  events  are 
required . 

The  Laplace  transform  of  the  output  inter-arrival  time  random 

variable  is  derived  as  follows:  Beginning  at  a  time  when  an  output 

event  has  Just  occurred,  let  Xq,  X2»’  . .  denote  the  time  between 

tii 

further  successive  input  events.  Assuming  that  the  N  event  produces 

the  first  count,  we  have  >  t,  Xg  >  t, ...jX^  >  t,  but  Xjj  <  t 

(note  that  there  are  no  conditions  on  Xq) .  Let  t(n)  =  xo+xi+*”+xn 

denote  the  total  time  between  output  events.  Then  the  conditional 

Laplace  transform  of  T(N),  given  N,  is 

N-l 

-s  X 

(3.2.1)  E[e'sT(N)|N]  =  E[e  °]  E[e  i=1|Xi  >  t]  E[e  n|x^  tl  . 


Noting  that  the  conditional  probability  that  two  input  events  are  X 
units  apart,  given  that  they  are  not  more  than  t  units  apart,  has 
density 


(5.2.2) 


Xe 


-XX 


1-e 


■Xt 
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we  find  the  conditional  Laplace  transform  (i.e.,  the  transform  of  this 
conditioned  random  variable)  is 


■«  „ 


(3.2.3)  E(e-BX|X  <  T)  -  f  '  e"eX  IX  *  - - — - 

Jo  1-e  (s+X)(l-e  ) 


Similarly,  if  we  know  that  two  input  events  differ  In  time  by  at  least 
r  units,  we  obtain  the  conditional  Laplace  transform 


(3.2.4) 


E[e'sX[X  >  t] 


“  -sX  Xe'U  Xe‘(s+X^ 

e  -xr  “  = , 


(s+X)e’ 


From  these  results  we  obtain  the  conditional  Laplace  transform  of 
T(N)  given  N 

m  “  SX  y  N-l  Y 

E[e”sT|N]  »  E[e  °](E[e'SX|X  >  t])  E[e'sX|x  <  t] 


(3-2.5) 


=  _K_  Xe-(s+X)T]N-l  X(l-e-(s+X)T 
X+s  (s+X)e"Xr  (s+X)(l-e’XT) 

,  N+l  -(N-l) (s+X)t  -(s+X)t 

/  a.  \  e _  x-e _ 

”  vX.+s; 


-(N-l)Xt 


1-e 


To  obtain  the  unconditional  transform  of  T(N)  we  multiply  by  the  pro- 

bability  that  the  first  output  event  will  occur  at  the  N  input 
-Vr  -X/r  N-l 

event,  namely,  (l-e  )  (e  )  and  sum  over  N  obtaining 


(3.2.6) 

This  yields 


E[e'sT]  =  l  E[e-Si  |N]  (l-e~A,T)(e"A’T)J; 
N=1 


-sTi 


■  Xt \ /  —  Xt \ N- 1 


(3.2.7)  E[e"sT]  =  (^$^)C  (1-e 


X  .2 


-(X+s)t 


)  2  f 

N-l 


Xe 


(s+X)t  N-l 


s+X 


-] 
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and  after  sunning  the  infinite  series 


(3.2.8) 


x.  2  ,  _-(X+s)T 

<P0(b)  =  (rrr)  - 

®  ,  _L. 

x  "  X,+s 


2  2 

In  this  formula,  the  term  X  /(X+e)  represents  the  transform  of  the 
waiting  time  for  the  second  input  event  of  the  Poisson  stream.  As 
T  -» oo  the  second  input  event  can  be  taken  as  producing  an  output 
event,  which  is  what  we  would  expect  of  an  infinite  memory  device. 

Proceeding  in  a  fashion  similar  to  that  used  for  the  Type  II 
counters,  we  obtain  after  some  manipulation, 


(3-2.9) 


<p '(s')  ir 


X(l-e 


-(X+SJTv 


(3.2.10) 


h  -  i  + 1] 


+  ^+s  r  1  ST  e~(^+s)T 

X  X(l-e"(B+X,T  "  X  (l.e-(Us)T)‘ 


d  r  1  1  If 
— o  Itrrn  J  ■  t  l- 


2  L?7^7J  =  X  \(l.e-(X+s)t. 


(5.2.11) 


+  Ja  ,  T  , 

x  y  (1_e-('T+s)tj2J 


so  that 


+  s  (  ...  3 


(3.2.12) 


1  +  1  2-e 

=  1  X(l-e’XT)  =  X(l-e-XT) 
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(3.2.13) 


_2  r  2-e  ,2  2 

0  X(l-e  Xx)  X2(l-e“XT) 

-Xt 

.  t  e 

+  I  * 

(1-e  ) 


(3.2.110 


ECM^T)}  ~  ^ 


(3.2.15)  varfM^T)}  ~  ^  [XTe'XT-2(l-e‘Xx)  +  (2-e‘Xx)2]  . 

Let  us  now  differentiate  E{M^(T) )  with  respect  to  X 


^  E(MX(T))  -  T(l-e"XT)(2-e’X’T)‘1 


(3.2.16) 


+  TXTe"XT(2-e”Xx)”X-TXTe’’Xx(l-e”XX) (2-e”Xx)-2  . 


Hence  the  efficiency  of  the  coincidence  counter  is 


C  =  |(T(l-e"X'r)(2-e"XT)"1  +  TXTe”xt(2-e'Xx)’1 


(3.2.17) 


.TXIe-XT(l-e-XT)(2-e-XT)*2)2  x 

f  (2-e'Xx)?  , 

XT(l-e"Xx)[XTe"Xx-2(l-e"^x)+(2-eAx)2] 


(3.2.18)  i 


(l-e"Xx)(2-e~Xx)"1+Xte’Xx(2-e~Xx)~1-XT(l-e~Xx)(2-e“XT)~2e~Xx]2 

(l-e‘Xx)2(2-e_XT)(TXe‘Xx-2(l-e'XT)+(2-e‘XT)2] 


Such  a  coincidence  counter  resembles  a  Type  I  counter  In  that  the 
expected  number  of  counts  E(M^(T)}  is  monotonically  Increasing  In  X. 

It  differs  from  this  type  of  counter  in  that  the  expected  number  of 
counts  is  asymptotically  proportional  to  the  number  of  input  events  as 
t  becomes  large.  Ujis  is  also  reflected  in  its  quantum  efficiency 
uhich,  although  more  complex  in  form,  tends  to  increase  from  zero  for 
X  »  0  to  100  percent  as  X  tends  to  infinity. 

The  response  function,  r(X),  of  the  coincidence  counter  is 
clearly  monotonic  since  formula  (3*2.16)  which  is  its  derivative  has 
no  zeros  in  (0,  oo).  Also,  as  is  apparent  from  the  expression  for 
£  ,  the  efficiency  of  the  counter  in  estimating  the  parameter  of  a 
Poisson  stream  is  independent  of  T. 

In  actual  practice  a  coincidence  counter  could  be  expected  also 
to  have  a  short  dead  time  after  registered  events,  so  that  a  cross  be¬ 
tween  a  Type  I  and  a  pure  coincidence  counter  would  have  to  be  considered. 

It  Is  also  interesting  to  speculate  as  to  the  suitability  of  such 
a  counter  for  separating  a  periodic  signal  from  contaminating  Poisson 
sources.  If,  for  example,  a  coincidence  counter  with  variable  time 
constant,  t,  could  be  constructed  and  a  periodic  pulse  with  repetition 
rate  Xq  were  Incident  as  well  as  a  Poisson  source  with  mean  rate  X, 
we  would  have  that 


(3.2.19) 


E(M(T))  »  XT 


for  t  a  little  less  than  and  X  «  X  (since,  in  combination 

o  0 

with  an  adjacent  periodic  pulse  every  random  pulse  produces  a  count) 
but 
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(5.2.20) 


*(M(T)) 


(X+^)T 

- 2 

for  t  >  y-  ,  so  that  a  discontinuity  in  the  plot  of  E(M^(T) )  as  a 
o 

function  of  -r  would  indicate  the  presence  of  the  periodic  pulse 
train  if  present.  (For  t  «  Xq  we  would  expect  that  periodicity 
would  be  irrelevant  so  that  in  that  case 


B(M(T)) 


-(X+X  )t 

(X+X  )  E(l-e 


2-e 


- Tzvrrr 


as  in  (3.2.14)).  , 

Coincidence  counters  have  been  constructed  and  used  (15)  to  study 
radiation  which  does  not  have  Boisson  character.  Here  again,  the  plot 
of  E{M^(T) }  for  fixed  mean  input  rate  X  and  fixed  observation  time 
T>  but  variable  time  constant  T  should  be  useful  as  a  comparison  with 
one  obtained  from  a  known  source.  One  can  use  it  to  infer  the  existence 
of  non-randomness  in  the  unknown  source  as  well  as  to  estimate  the 
parameters  of  interest. 


3.  'Pxe  R-Fold  Coincidence  Counter 

One  natural  extension  of  the  simple  coincidence  counter  is  to  per¬ 
mit  input  events  arriving  within  specified  time  intervals  to  change 
state  levels  of  a  counter,  and  allow  the  registering  of  an  output  event 
only  after  a  terminal  state,  say  after  the  R^*1  input  event,  is  reached. 
We  term  this  the  R-fold  coincidence  counter.  For  R  «  2,  this  reduces 
to  the  previously  described  simple  coincidence  counter  which  has  two 
states  S  and  S. .  Ihen  an  input  event  changes  the  counter  from  S 
to  8^,  and  if  another  event  occurs  before  time  t,  an  output  count  is 
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registered  and  the  counter  immediately  reverts  to  SQ.  The  counter  also 
lmaedlately  reverts  to  SQ  at  time  t  without  registering  If  the  Input 
does  not  arrive  within  the  time  interval. 

Briefly  we  can  describe  the  R-fold  coincidence  counter  as  follows: 

A  counter  which  immediately  registers  an 
output  event  whenever  an  input  event  is 
preceded  by  (R-l)  non-registered  input 
events  none  of  which  are  spaced  more  than 
x  units  apart  and  r  is  a  time  interval 
of  random  length  whose  (R-l)  values  are 
independently  and  identically  distributed. 

Here  we  shall  derive  the  waiting  time  distribution  (or,  rather,  its 
Laplace  transform)  for  the  first  time  occurrence  of  the  critical  level  R 
in  the  case  where  the  waiting  time  between  upward  jumps  of  unit  magnitude 
has  negative  exponential  distribution,  F(t)  =  1-e-^.  This  is  equiva¬ 
lent  to  the  assumption  that  input  impulses  constitute  a  Poisson  process. 
We  also  asstime  that  all  the  levels  Sy  J=l,2, . ..,R-1  can  decay  only 
to  Sq  and  that  such  decay  takes  place  at  the  end  of  an  interval  of 
duration  following  the  latest  incident  impulses.  We  asstime  that 
each  tj  is  itself  a  random  variable  with  absolutely  continuous  density, 
and  with  distribution  function  G(r)  which  is  independent  of  Sj  and 
the  times  of  occurrence  of  the  states  S  .  The  special  case  R  *  2; 

G(t)  =  0  for  t  <  tq,  G(r)  =  1  for  t  >  x ,  was  treated  by  another 
method  in  the  previous  section. 

We  shall  not  attempt  to  evaluate  the  efficiency  of  such  a  general 
coincidence  counter  explicitly  since  the  mathematical  manipulations  be¬ 
come  too  tedious  and  unmanageable.  Ad  hoc  cases  of  special  Interest 
can  be  handled  by  computers.  Our  concern  in  this  section  is  primarily 
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to  arrive  at  the  formula  for  the  inter-arrival  time  distribution  from 
which,  by  the  methods  covered  in  earlier  sections,  the  efficiency  can 
be  derived  for  those  specific  cases  that  may  be  of  interest. 

The  process  with  which  we  are  concerned  is  non-Markovlan  in  the 

sense  that  merely  knowing  which  state  8j,  J«1,2,...,R-1,  the  process 

is  in  at  a  given  time  t  is  not  sufficient  to  calculate  the  probable 

behavior  of  the  process  beyond  this  time.  In  addition,  we  need  to 

know  how  long  ago  the  last  impulse  occurred.  The  only  case  in  which 

such  knowledge  is  not  needed  is  when  at  time  t  the  process  is 

discovered  in  state  S  . 

o 

We  shall  write  Pj(t,r)  for  the  probability  density  that  Jointly: 
a)  at  time  t  the  process  is  in  state  Sj,  b)  that  it  reached  this 
state  t  time  units  ago,  that  is,  at  t  -  t;  c)  that  the  state  at 
t  =  0  is  Sq  ,  and  that  previous  to  time  t  the  state  R  was  not 
reached.  Now  the  probability  that  a  decay  to  SQ  takes  place  during 
the  small  interval  (t,t  +  6)  after  occurrence  of  the  last  impulse  is 
8  G'(t),  and  the  probability  that  no  decay  takes  place  during  the  time 
interval  (t-T,t)  is  1-G(t).  Hence  the  conditional  probability  of 
decay  during  the  interval  (t,t  +  8),  given  that  no  decay  has  taken 
place  during  the  elapsed  time  t  since  the  last  impulse,  label  it 
8.p(x),  is  given  by 

(3.3.1)  8.p(x)  =  6‘G>(T^  . 

1  -  G(t) 

Later  on,  it  will  be  useful  to  consider  the  solution  of  this  differential 
hence  we  note  here  that  from  the  above  one  can  obtain 

3* 


equation, 


(3.3.2) 


- /  n(x)dx 
G(t)  -  l-e  Jo 

The  function  u(f )  is  sometimes  called  the  hazard  function,  and  has 
been  discussed,  for  example,  in  connection  with  telephone  call  demands 
and  other  queuing  problems.  In  the  particular  case  where  G(t)  is  a 
negative  exponential  distribution,  it  is  seen  that  u("0  is  constant, 
so  that  the  conditional  probability  of  decay  after  t  time  units  from 
the  last  impulse  have  elapsed  is  independent  of  t.  It  is  in  that 
sense  that  the  negative  exponential  distribution  is  said  to  have  no 
memory. 

We  shall  derive  here  a  set  of  differential  equations  from  which  the 
Laplace  transform  inter-arrival  time  distribution  of  output  counts  can 
be  obtained.  First  let  us  obtain  the  probability  that  the  state  of  the 
process  will  be  Sj  (j  >  0)  at  time  (t+ft)  and  that  t+6  time  units 
have  elapsed  since  the  last  impulse,  given  that  at  time  t  the  time 
elapsed  since  the  last  impulse  was  t  and  write  it  as  a  function  of 
Pj  (t, t ) •  We  assume  here  that,  for  small  6,  the  probability  of  an 
impulse  during  a  time  interval  8  is  X8  and  the  probability  of  no 
impulse  is  1-X8}  and  that  other  cases  have  such  small  probability 
that  they  can  be  neglected.  These  are  the  usual  assumptions  linked 
to  the  fact  that  the  input  impulses  form  a  Poisson  stream.  In  terms 
of  Pj(t,t)  and  the  hazard  function  u(t),  we  can  write 

(3.3.3)  Pj (t+8,T+6)  -  (1-X8)(1-5h(t))pj (t,r) 


35 


or,  neglecting  terms  in  &  slnoe  8  is  assumed  small  we  get 


(3.3.U) 


Pj(t+8,T+8)  -  [1-X8-6  u(T)3Pj(t,T)  . 


If  during  the  interval  (t,t+6)  an  impulse  occurs,  we  have,  for 

the  probability  that  the  state  is  S.  at  t+6  and  that  the  time  t' 

0 

since  the  last  impulse  was  lees  than  8  [denote  this  by  p^ (t+8,t'<  8)] 

r  co 

(3.3.5)  Pj(t+8,t'<8)  -  X8  J  (l-6li('r))pJ_1(t,T)dT 

where  the  integrand  on  the  right  represents  the  probability  that  the 
state  was  Sj  ^  at  time  t  ,  had  been  attained  t  time  units  earlier; 
and  did  not  spontaneously  decay  in  the  time  span  (t,t+B)  following 
the  moment  it  was  attained.  Since  all  possible  times,  T,  of  attain¬ 
ment  of  the  state  S.  T  must  be  considered,  the  probability  of  no 

f  q° 

decay  from  state  Sj_1  during  (t,t+&)  is  J  (l-Bu(T))pj_1(t,T)dT. 
Again  neglecting  terms  in  82  this  becomes 

roo 

(3.3.6)  Pj(t+8,T'<  8)  =  18  J  pj^Ct^JdT  . 
r°o 

Note  that  J  p^  1(t,T)d-r  is  the  probability  that  the  state  is 

S  at  time  t  regardless  of  how  long  ago  the  state  was  attained. 

For  use  later  on  we  note  that  we  can  approximate  Pj(t,8,T'<  8)  by 

B.p.(t,6,T'  =  0)  since  the  time  interval  8  is  small,  so  that  the 
J 

probability  of  an  impulse  at  some  instance  t'  in  (0,8)  is  nearly 
the  same  for  all  t'  in  this  interval  [see  (3.3.20)]. 


Itaeae  considerations  hold  for  j-1,2, . . .,R-1.  For  state  Sq,  we 

have 


R  X  oo 

(3.3,7)  p  ft+6)  -  (1-18 )p(t)  +  j  (1-15)  f  6  |i(T)p,(t,t)dT 
°  °  J-l  Jo  J 


l.e.,  the  probability  of  being  in  state  SQ  at  time  t+8  is  the  pro¬ 
bability  of  being  in  that  state  at  time  t  and  having  no  incident 
events  in  (t,t+B)  plus  the  probability  that  a  decay  occurs  in  the 
period  (t,t+8),  when  the  process  is  in  state  ( j=l,2, . . ,,R-l)  at 

time  t  and  no  incident  events  occur  during  (t,t+B).  Note  that  the 
integral  under  the  summation  represents  the  probability  of  a  decay  in 
the  period  (t,t+8)  when  the  process  is  in  Btate  S  (j=l,2, . . .,R-l) 
at  time  t  and  this  probability  is  independent  of  how  long  ago  the 
state  Sj  was  attained  previous  to  time  t.  Similarly  for 
Po(t+5,T'<  8)  we  have  by  (3.3.5) 

r00 

(3.3.8)  Po(t+B,T'<5)  =  18  J  (l-ButT^p^^OdT  . 

In  the  period  (t,t+8)  a  transition  from  state  SR  ^  to  state 
may  occur  thus  causing  the  process  to  revert  immediately  to  state 
Sq  and  register  an  output  count,  or  we  may  write 

R**l  oo 

(3.3.9)  P„(t+8)  =  (l-X8)p  (t)  +  i  [  MOp.(t,t)dT 

J-lJo  J 


and 


(3.3.10) 


P0(t+5,T'< 


Po(t+8,8')d8' 


r00 

IB/  pR_-(t,T)dT 

Jo 
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2 

if  terms  In  &  are  neglected.  If  In  (3.3>*0  we  transpose  Pj(t,x), 
then  divide  each  side  of  the  equation  by  &  and  pass  to  the  limit, 
we  obtain 


lim 

&-»0 


p  (t+6,T+6)-p.(t,t) 


lim 

&->0 


or 


dp  dp 

(3.3.11)  *  -(X.  +  l^(  T )  )  Pj(t,x)  . 


lhis  is  so  for  J.*l,2, . . .,R-1  and  similarly  from  (3.3.9)  we  get  for 
J  -  0 


(3-3.12) 


dp  R-l  r  00 

fir  *  -X.  P0(t)  +  2,^  J  |l(T)Pj(t,T)dT  . 


These  are  the  differential  equations  whose  solutions  will  provide  the 
inter-arrival  time  distribution  of  the  output  counts.  First  we  will 
obtain  an  expression  for  the  probability  of  attaining  state  R  at 

rT 

some  point  in  [0,T],  I  Pp(t')dt'.  From  this  we  find  the  probability 
density  for  the  waiting  time  until  state  R  is  reached.  This  is 

rT 

obtained  by  differentiating  I  pR(t'  )dt'  with  respect  to  T,  since 

Jo 

this  can  be  interpreted  as  the  probability  that  the  waiting  time  for 
reaching  state  R  will  not  exceed  T.  The  remainder  of  this  discussion 
is  therefore  concerned  with  obtaining  this  distribution  by  developing 
its  Laplace  transform. 

Equation  (3.3.11),  treated  as  a  differential  equation  in  the 
variable  t,  can  be  solved  by  taking  the  Laplace  transform  with 
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respect  to  the  variable  t.  Label  (b,t)  the  Laplace  transform  of 


Pj(t,-r)  or 


(3.3.13) 


■  1 


oo 


«"8^  Pj(t>x)dt  . 


Nov  take  the  Laplace  transform  of  each  side  of  (3.3.11)  and  (3.3.12)  and 
note  that 


(3.3.14) 


r  <ju  dp 

/  e  d r  dt  =  -  po(0)  +  8  n0(8) 

Jo 


with  PQ(0)  =  1,  while 


(3.3.15)  f  e 
Jo 


00  -St  * 


3*4  dt  =  -  p j ( 0, t )  +  s  «j(b,t)  =  s  *j(s,t) 


with  Pq(0,t)  =  0  for  j  /  0. 

Hence  equation  (3. 3. 11)  becomes 


(3.3.16)  s  «j(s,t)  +  ^  «j(s,t)  =  -  [x  +  u(t)]itj(s,t) 


and  thus  by  rewriting  we  get 

(3.3.17) 


«j(s,t)  =  nj (s,0)  e" (B+^)  e 


-f 

T  JO 


n(x)dx 


or,  using  (3.3.2) 


(3.3.18)  *j(s,t)  =  itj (s,0)  e"(8+X)T[l-G(x)]  . 


Since  we  are  concerned  only  with  the  marginal  distribution  of 
Pj(t,t),  we  integrate  over  t  to  obtain  [setting  1-G(t)  »Q(t))  and 
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and  using  (3.3.l)>  this  becomes 


(3.3.25)  J  n^(B,x)n(r)di  =  Jtj(e,0)  J  e‘^8+X^T  G'(t)dT 


«j (s,0)  G'*  (s+X)  . 


Finally,  we  take  the  Laplace  transform  of  (3.3*12). 

R-l  f 00 

(3.3.26)  -1  +  8  it  (s)  =  -X  it  (s)  +  £  / 

j=l  Jo 

and  substituting  the  result  (3*3.25)> 


Kj(s,t)u(t)  dr 


R-l  _ 

(3.3*27)  -1  +  an  (s)  =  -X  Jt  (s)  +  J  *^(8*0)  G'  (s+X) 

j=l  J 


=  -Xrt 


R-1  ,  .  . 

,(s)  +  «Q(s)  O'  (s+X)  X  X1  (y(8+^))J 


On  summing  we  obtain  the  result 


(3.3.28)  -1  +  s  n  (s)  =  -X  n  (s)  +  G'  (s+X)  X  n  (s) 

00  °  1  1-X3*(s+X) 


i-(x  a*(s+x)) 


R-n 


or. 


(3.3.29) 


s  +  X-X  G1’ 


L  1-X^»(s+X) 


If  we  let  H(t)  denote  the  probability  that  a  count  (that  is, 
state  SR)  is  attained  in  the  interval  (0,T),  we  have  from  the 
fact  that 


4l 


(3.3.30) 


rT 

H(T)  -  J  PpCt'Jdt' 


or,  on  different dating. 


(3.3.31) 


H'(T)  =  pR(T) 


that 


(3.3.32) 


-8tp 

E[e  R)  -  (H'(T))*  =  nR(s) 


is  the  desired  Laplace  transform  of  the  waiting  time  distribution 
towards  which  we  have  been  working.  For  state  SR  ,  we  also  have 


(3-3.33) 


jtr(s)  =  «R(s,0) 


since  the  system  immediately  returns  to  state  SQ.  Hence,  using 
(3.3.21), 


rOD 

(3. 3. 3*0  *R(s)  m  x  J  5tR-i^B»T)dT 

and  substituting  from  (3.3.19), 


(3.3.35) 


nR(s)  -  \  *^(8,0)  !*(s+X.) 


or.  Iterating  this  substitution  as  in  (5.3.25); 


(3.3.36) 


JtgU)  -  {^*(s+X))R_1  no(s) 


k2 
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CHAPTER  IV 


QUANTUM  EFFICIENCY  OF  KOTOGRAMIC  DETECTORS 

1 .  On  the  Problem  of  Specifying  a  Measure  of  Performance 

In  one  sense,  the  perfect  photographic  detector  is  one  in  which 
only  a  single  incident  photon  is  required  to  make  a  grain  developable. 
Consequently,  if  an  average  of  XT  photons  is  incident  per  grain 
during  the  time  of  exposure,  T,  and  if  the  Poisson  distribution  ap¬ 
plied,  the  probability  of  a  grain  being  made  developable  is  the  sum 
of  the  probabilities  of  1,  2,  3,..*  photons  incident  on  a  grain  i.e., 

(4.1.1)  H  [0,T]  =  H(XT)  =  £  e"XT  =  l-e'XT  . 

k  j=l 

An  "imperfect"  emulsion,  however,  can  be  constructed  in  at  least  three 
ways. 

In  the  first  place,  the  grains  may  not  be  tightly  packed.  In  this 
case,  even  though  the  flux  in  photons  per  grain  is  XT  and  a  fraction 
H(XT)  of  the  grains  in  a  tightly  packed  emulsion  could  be  made  develop¬ 
able,  only  a  fraction  Q-H(XT)  are  actually  made  developable,  where 
Q  is  the  ratio  of  the  number  of  grains  present  to  the  number  that 
could  be  packed  into  the  same  surface  area.  If  only  a  single  photon 
is  required  for  developability,  the  fractional  number  of  grains  made 
developable  in  such  a  case  is 

(4.1.2)  Q'H(XT)  -  Q  £  e‘XT  -  Q(l-e‘XT) 

In  the  second  place,  the  grains  may  not  be  good  absorbers,  so 
that  the  effective  flux  per  grain  is  not  XT  but,  say,  QXT  ,  where 
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Q  now  denotes  the  absorption  of  a  grain.  This  is  equivalent  to  as sum- 
ing  perfect  absorption  and  placing  a  filter  of  transmission  Q  over 
the  emulsion.  Then  the  number  of  grains  made  developable  (assuming  one 
photon  per  grain  required  for  developability)  would  be 


(4.1.5)  H  (\T)  *  l  .-««  . 

W  J«1 

However,  it  is  also  possible  to  think  of  the  photographic  detector 
as  being  imperfect  in  the  sense  that  a  grain  might  require  R  =  l/Q 
photons  to  become  developable,  (in  this  case  we  can  think  of  Q  as  a 
measure  of  the  inertia  of  a  grain.)  If  each  grain  were  so  constituted, 
the  fraction  made  developable  would  be  the  sum  of  the  probabilities  of 
R,  R  +  1,  R  +  2, ...  photons  incident  on  a  grain,  i.e.. 


(4.1.4) 


Hr(XT) 


=  l 

J=R 


a-XT  (XT) 


J! 


Formulas  (4.1.2)  -  (4.1.4)  indicate  the  three  quite  different 
ways  in  which  a  "fractional  utilization  factor,  Q"  could  be  introduced. 
In  general,  an  emulsion  will  be  imperfect  as  a  result  of  all  these 
reasons.  For  a  combination  of  the  last  two  reasons  we  could  write 


(4.1.5)  Wxt)  .  jB  .-«*  sagt 

or,  even  more  generally, 


(4.1.6) 


\,Wm 


I  0,  I  e 

i-i  ^ 


-(QW) 


(QXT)J 

Jl 
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In  ease  the  emulsion  consists  of  a  variety  of  grains  and  the  fraction 
of  the  population  requires  photons  for  developability.  In 
formula  (4.1.6),  represents  the  fractional  area  covered  by  grains 
of  inertia  when  all  the  grains  absorb  a  fraction  Q  of  the  In¬ 
cident  energy.  (Holes  are  considered  as  grains  with  infinity  inertia.) 

Aside  from  the  fact  that  "fractional  utilization"  factor  is  an 
ambiguous  notion  as  shown  above,  there  are  other  reasons  why  the  old 
definition  of  quantum  efficiency  as  a  ratio  of  number  of  output  to 
input  events  (which  is  linked  to  the  notion  of  fractional  utilization) 
is  not  a  good  measure  of  performance. 

If  every  input  event  gave  rise  to  a  single  count  we  should  like 
to  say  that  the  detector  has  a  quantum  efficiency  of  unity.  Even  if 
every  input  event  gave  rise  to  several  counts  we  should  still  wish  to 
speak  of  the  device  as  having  unit  quantum  efficiency,  and  we  should 
consider  each  aggregate  of  counts  produced  by  one  input  event  as  a 
single  Output  event.  Suppose,  however  that  on  the  average,  only  every 
other  input  event  is  effective  in  producing  output  events  and  that  for 
each  output  two  counts  are  registered.  The  ratio  of  output  to  input 
counts  ia  one  and  yet  a  detector  in  which  every  input  event  produces 
an  output  count  may  be  viewed  as  superior.  It  follows  that  a  different 
approach  to  a  unique  and  meaningful  definition  of  quantum  efficiency 
is  called  for— one  which  cannot  only  be  computed  theoretically  but 
can  also  be  measured  experimentally.  Ibis  is  accomplished  by  applying 
to  photographic  detectors  the  concepts  developed  in  Chapter  I. 
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*  Bi£  SSSSSB  Bfflcig»cy  of  Simple  "Photographic"  X-Detectors 

We  shall  consider  a  (simple)  photographic  detector  as, an  ensemble 
of  go-no-go  detectors,  each  of  which  is  capable  of  responding  just 

XV 

once  immediately  following  the  arrival  of  the  R  event.  Consequently, 
further  arrivals  at  a  detector  having  already  received  R  "hits"  are 
wasted,  and  cannot  be  transferred  to  another  detector  or  registered 
in  the  output.  The  photographic  detector  is  unique  in  that  it  can 
operate  simultaneously  on  a  large  set  of  Poisson  processes,  namely, 
on  all  those  sources  in  space  which  are  imaged  on  the  face  of  the 
photographic  plate.  Also,  because  of  imperfect  imaging  and  scatter 
within  the  emulsion,  a  Poisson  point  source  is  imaged,  not  on  one, 
but  on  a  set  S  of  detectors  (or  photographic  grains).  Thus,  in 
comparing  the  intensity  of  two  point  sources,  the  output  from  2S 
detectors  must  be  compared.  An  ideal  non-photographic  ensemble  of 
S  detectors  exposed  to  a  Poisson  source  of  intensity  X  for  a  time 

A  X 

T  would  provide  an  estimate  X,  of  X  with  a  variance  of  ^  ,  assum¬ 
ing  the  set  of  S  detectors  to  operate  independently  of  each  other  on 
the  same  Poisson  source. 

If  we  suppose  that  each  photographic  detector  requires  exactly  R 
hits*  to  become  developable  and  retains  its  developability  indefinitely, 
then  we  can  calculate  that  the  mean  number  of  detectors  responding  to 
an  intensity  X  (per  grain  per  unit  time)  after  an  exposure  time  T 


*7 


will  be 


(4.2.1) 


S.H(XT)  -  S[  £  e"XT 

A  3* 


Now  If  H  Is  the  probability  that  a  detector  will  respond,  and 
if  S  independent  experiments  are  performed,  the  probability  that 
exactly  k  responses  are  obtained  is 

(4.2.2)  p(k)  -  (®)  (l-H)S‘k 

and  hence  the  variance  of  the  number  of  responses  is 

(4.2.3)  =  S*H(XT) [I-H(XT) ]  . 

Since  (4.2.1)  yields  the  response  function  r(XT)  of  the  photo¬ 
graphic  detector  (known  in  photographic  theory  as  the  characteristic 
curve) and  (4.2.3)  the  variance  of  the  output,  we  can  use  (1.3.5)  to 
obtain  the  variance  of  the  asymptotically  unbiased  estimate  of  XT  , 
r_1( XT),  so  that 

(4.2.4)  var{x|M(T))  -  [^Ly.]  ‘2  . 

As  indicated  earlier,  the  ratio  of  this  variance  to  that  of  the  ideal 
detector  is  a  measure  of  the  efficiency  of  the  detecting  process. 

To  illustrate  the  calculation  of  quantum  efficiency  with  a  simple 
example,  suppose  that  we  have  a  photographic  detector  requiring  one 
photon  hit  to  make  each  grain  developable.  If  the  total  radiation  on 
one  grain  during  the  time  of  exposure  has  average  XT,  the  fractional 
number  of  grains  which  will  become  developable  is 


(4.2.5) 


H(XT)  .  I  . 

j-1 

It  there  are  a  total  of  S  grains  within  a  given  area,  we  expect 
there  to  he 

(4.2.6)  S  [H(XT)  ]  =  S(l-e"XT) 

developed  grains  after  exposure  and  processing.  Due  to  random  varia¬ 
tions  in  the  incident  events,  there  will  be  fluctuations  in  the  number 
of  developed  grains  from  one  area  of  S  grains  to  the  next.  Hie 
standard  deviation  of  the  number  of  developed  grains  is  obtained  from 

/~  \m  ™  \ m  \ffl 

S(l-e"  )e"  (using  1-e”  as 
the  probability  of  success).  The  estimate  of  X  is 

(4.2.7)  X  =  H-1(XT)  =  -  |  log(l-H(XT)) 

and  the  variance  of  the  estimate  is  by  (1.5-5) 

var(X)  =  var(S[H(XT)])  -] 

(4.2.8) 

Se-Xtd-e-XT) 

=  S2lV2XT 

whereas  the  variance  of  X  in  the  case  of  a  set  of  S  ideal  non¬ 
photographic  detectors  is 

(4.2.9)  var(X)  =  ~  . 

Thus  the  quantum  efficiency  is 
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(4.2.10) 


c 


5  (XT, R-l)  -  ££ 


ST2  e’M 
(l-e'^T)e-^ 


XTe' 


•XT 


1-e 


Tt 


> 


and  this  depends  on  the  total  expected  number  of  input  events  XT 
rather  than  on  the  input  rate  X.  We  also  see  that  d  is  a  monotonically 
decreasing  function  of  XT,  and  that  the  one  photon  photographic  de¬ 
tector  has  100$  efficiency  only  at  XT  -  0.  For  given  X,  this  de¬ 
tector  becomes  increasingly  inefficient  as  T  increases.  The  photo¬ 
graphic  type  of  detector  is  a  saturating  device  which  has  infinite 
memory. 

In  general  the  photographic  detector  requiring  R  photons  for 
developabillty  has  its  peak  quantum  efficiency  at  about  XT  =  R-l. 

This  follows  from  the  fact  that  dfl/d(XT)  is  a  Gamma  function  in  XT 
with  parameter  R-l,  and  hence  its  mode  is  at  R-l.  It  follows  that 
if  the  background  radiation  is  known  to  have  mean  value  XT,  the  ideal 
photographic  detector  is  the  one  vhich  requires  R  =  XT  +  1  photons 
for  developability  rather  than  one  photon.  At  the  value  of  XT  *  1,  the 
efficiency  of  the  R-photon  detector  haB  the  form 


(4.2.11) 


R-l 

H(l-H) 


-(R-l),  vR-1  2 


Using  Stirling's  formula  for  nl 

n  +  §  -n 

(4.2.12)  nl  ~  n  d  e  n 


we  obtain 


a 
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(4.2.13) 


The  maximum  value  of  H(l-H)  Is  .25,  so  that  a  lower  bound  for  the 
quantum  efficiency  peak  value  Is 

Independent  of  the  number  of  photons,  R,  required  for  developability. 
Ibis  is  perhaps  the  most  surprising  result  concerning  this  class  of 
radiation  detectors. 

When  the  detection  problem  is  one  of  detecting  the  smallest  pos¬ 
sible  incremental  signal  for  a  given  radiation  background,  the  corres¬ 
ponding  R-photon  emulsion  must  be  considered  as  the  ideal  photographic 
detector;  the  efficiency  of  an  actual  eraulBion  relative  to  this  detector 
is  obtained  by  dividing  the  maximum  quantum  efficiency  of  the  actual 
emulsion  by  the  quantum  efficiency  of  the  appropriate  R-photon  photo¬ 
graphic  detector.  This  yields  a  somewhat  higher  value  of  efficiency 
than  appears  in  the  literature. 

However,  even  the  higher  value  of  efficiency  obtained  ip  this 
manner  is  not  a  really  fair  index  of  improvability  for  actual  photo¬ 
graphic  materials.  The  reason  for  this  is  that  the  assumption  of  a 
given  background  radiation  is  unrealistic.  We  do  not  know  the  strength 
of  the  background  radiation.  This  means  that  the  photographic  detector 
must  have  a  certain  latitude.  The  latitude  we  require  is  a  measure  of 
our  a  priori  ignorance  as  to  the  strength  of  the  radiation  background. 
All  we  can  do  in  the  face  of  this  ignorance  is  to  make  the  quantum  ef¬ 
ficiency  as  uniform  as  possible  for  the  range  of  background  radiation 
strengths  that  we  expect  to  encounter.  If,  for  example,  we  expect  a 
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range  of  1  to  10  photons  per  grain  area  for  the  background  radiation, 
we  might  take  for  our  ideal  photographic  detector  a  composite  of 
R  -  1, ...,8.  Such  a  composite  detector  made  up,  say  of  an  equal  number 
of  grains  of  each  type  has  a  more  uniform,  but  lower,  quantum  efficiency 
than  a  one-photon  detector  would  have.  If  a  still  greater  latitude  is 
required,  say  for  a  range  of  XT's  from  1  to  100,  a  uniform  quantum 
efficiency  over  this  range  can  be  achieved  only  by  using  a  still  greater 
range  of  R's. 

If  we  want  to  determine  the  optimum  response  curve  over  a  given 
range,  (XQT,  X^T)  for  uniform  quantum  efficiency  throughout  this 
range,  we  proceed  as  follows:  Since 

*  h(i-h) 

and  since  we  require 

(4.2.16)  (f  =  constant  XqT  <  XT  <  X.jT 
we  obtain 

1  11 

(4.2.17)  H'(XT)  (H(XT)[1-H(XT)])5/(XT)^  . 

This  leads  to 

11  1 

(4.2.18)  H(XT)  =■  sin2  (g?[(XT)5  -  (XqT)5]} 

1  1 

3  p  p 

with  XqT  <  XT  <  [n/2  £  +  (XqT)  ]  for  the  response  curve  having  the 

optimum  shape.  With  XqT  »  1  and  a  uniform  quantum  efficiency  of  1# 
n2 

a  value  X^T  “  -250  is  obtained.  The  range  is  thus  1  to  250. 
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Conversely,  if  the  range  la  restricted  to  the  values  X.qT  -  1,  X^T  ■  100, 
the  greatest  possible  uniform  quantum  efficiency  is  found  to  be  about 
2.5*. 

It  follovs  from  these  considerations  that  actual  emulsions,  which 
have,  at  present,  a  peak  quantum  efficiency  of  between  0.5*  and  1.0*  , 
differ  from  the  attainable  optimum  by  a  factor  which  is  certainly  less 
than  10.  The  improvability  of  photographic  materials,  given  the  con¬ 
straints  imposed  by  pictorial  photography,  is  thus  much  smaller  than 
the  factor  of  100  which  was  believed  possible,  and  no  further  dramatic 
breakthrough  in  photographic  "speed"  should  be  anticipated. 

3.  The  Photographic  Detector  with  Reciprocity  Failure 

As  in  the  previous  section,  we  postulate  that  a  photographic 
detector  consists  of  an  ensemble  of  go-no-go  detectors.  Once  one 
of  these  detectors  has  acted  (i.e.,  once  a  silver  halide  grain  has 
become  developable),  further  incident  hits  are  wasted.  Reciprocity 
failure  means  that  the  photographic  detector  responds  not  just  to  the 
total  number  of  photons,  XT,  incident  during  the  time  of  exposure, 
but  reacts  differently,  depending  on  whether  for  XT  *  constant  it  is 
the  time  of  exposure  or  the  strength  of  radiation  which  is  large. 

There  are  two  types  of  reciprocity  failure}  high  intensity  failure, 
which  can  be  attributed  to  a  type  of  dead  time  phenomenon  as  in  a 
Type  I  counter,  and  low  intensity  failure,  which  can  be  thought  of  as 
due  to  the  finite  memory  of  a  coincidence  type  of  counter  mechanism. 
Diverse  explanations  of  this  phenomenon  have  been  given  (9)  but  the 
crucial  experiments  to  determine  the  precise  mechanism  whereby  reciprocity 
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failure  is  produced  have  not  yet  been  attained.  Consequently  any  model 
whose  consequences  are  In  reasonable  accord  with  existing  experimental 
data  can  provide  a  step  forward. 

We  assume  here  that  for  photons  in  the  visible  region  there  is 
not  sufficient  energy  to  produce  a  developable  grain.  As  is  well 
known,  at  shorter  wave  lengths  there  is  enough  energy  in  a  particle 
to  trigger  one  or  more  photographic  grains,  and  for  such  particles 
the  problem  of  reciprocity  failure  does  not  arise. 

For  low  intensity  reciprocity  failure,  we  assume  that  a  photon  can 
produce  an  imbalance  in  a  silver  halide  crystal  which  can  persist  for 
a  time  t.  Another  photon  incident  during  this  time  will  cause  this 
imbalance  to  cease  by  producing  an  atomic  silver  speck.  Such  a  speck 
may,  or  may  not,  in  itself  be  developable.  If  it  is,  we  have  a  two- 
photon  photographic  detector  with  low  intensity  reciprocity  failure. 

If  t  Is  infinite,  this  reduces  to  the  type  of  photographic  detector 
discussed  in  Section  2.  It  may  be  necessary  to  have  a  larger  speck 
of  silver  to  produce  development  than  one  obtained  from  two  photons. 

If  we  assume  that  the  two  photon  speck  is  non-developable  but  one 
twice  as  large  is  developable,  then  two  further  photon  hits  within 
an  interval  T  are  required  to  produce  either  another  speck  or  to 
enlarge  the  one  already  formed.  Since  actual  emulsions  are  a  mixture 
of  grains  of  varying  sensitivity,  we  would  have  to  combine  various 
models  to  simulate  an  actual  photographic  material.  Here  we  shall 
content  ourselves  with  some  discussion  of  the  two  and  four  photon 
photographic  detectors,  as  even  these  present  considerable  difficulty. 
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To  arrive  at  the  response  curve  and  quantum  efficiency  of  a  two 
photon  photographic  detector  with  low  Intensity  reciprocity  failure, 
we  should  proceed  as  follows: 

Beginning  with  the  Laplace  transform  of  the  inter-arrival  time 
distribution  of  the  coincidence  counter  input  (3.2.8),  we  find  the 
inverse  Laplace  transform  and  integrate  this  transform  from  zero  to 
T.  The  result,  Fq(X,T,t),  indicates  the  probability  of  one  or  more 
coincidence  events  which  is  the  probability  of  a  grain  becoming  deve¬ 
lopable  during  the  time  of  exposure  T,  that  is,  the  probability  that 
at  least  one  silver  speck  is  formed  in  a  grain  irradiated  by  an  average 
of  XT  photons,  (and  one  such  speck  is  sufficient  for  developability) . 
If  there  are  S  photographic  grains  in  the  area  under  consideration, 
the  expected  number  of  grains  which  will  contain  a  developable  speck 
after  exposure  time  T  will  then  be  S.Fq(X,T,t).  A  plot  of  Fo(X,T,t) 
versus  X  or  T  indicates  the  average  fractional  number  of  grains 
which  become  developable  as  X  or  T  increases,  the  other  variable 
being  held  constant,  (in  photographic  technology,  when  studying  reci¬ 
procity  failure,  it  is  customary  to  hold  Fq(X,T,t)  (which  corresponds 
to  the  developed  optical  density) constant  and  plot  for  various  values 
of  X  or  T  the  value  of  XT  needed  to  produce  a  fixed  Fq(X,T, r). 

Once  we  have  Fq(X,T,  t)  we  can  obtain  the  variance  in  the  number 
of  developed  grains  from  the  binomial  distribution,  that  is,  the  vari¬ 
ance  in  the  number  of  developable  grains  can  be  expected  to  be  . 
SFo(X,T,t)*(1-Fq(X,T,t)).  Since  there  are  SXT  photons  expected 
during  the  time  of  exposure,  an  ideal  detector  would  estimate  X  with 
variance  X/ST  so  that  the  efficiency  of  the  photographic  detector  is 


(k.3.1)  £  -  ^  [SSglSl]2 


If  instead  of  two  photons  it  takes  four  photons  for  developability, 
we  must  convolve  F  with  itself  to  obtain  the  probability  of  at  least 
two  output  events  in  the  coincidence  counter. 

The  problem  of  carrying  out  this  procedure  arises  right  at  the 
start  in  trying  to  obtain  an  explicit  closed  form  for  the  inverse 
Laplace  transform.  To  find  the  inverse  transform  of  the  coincidence 
counter  inter-arrival  time 


(4.3.2) 


\  2  ,  -(X+s)t 

*0{e)  *  "  X  -(X-fs)T 

1  "  X.+8  ® 


one  should  take  a  contour  integral  over  the  left  half  plane.  (From 
the  fact  that  the  distribution  function  is  zero  for  t  <  0,  we  know 
that  all  the  poles  of  <p(s)  must  have  negative  real  part  (9).  Hence 
we  must  find  the  roots  of  s  +  X(l-e”^+8^T) . 

To  find  these  let 


(4.3.3) 


s  =  re  =  r  cos  6  +  ir  sin  9 


so  that 


(4.3.4) 


-T,s  --rr  cos  9  -rrisin  9 
e  =  e  e 


(4.3.5)  e"irT  008  6  =  cos(rT  cos  9)  -  i  sin(rr  cos  0) 


so  that 


(4.3.6) 


+  V[l-e-(8+X)T]  =  0 
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becomes 


(4.3.7)  I.  r  cos  9  +  X-Xe”TX,’Tr  008  0  cos  (it  cob  -  0 
and 

(4.3.8)  II.  r  sin  9  -  Xe”**  e”Tt  008  0  sin  (it  cos  9) l  -  0  . 

Squaring  and  adding  these,  we  obtain 

(4.3.9)  r  =  +  X[l-e  e  ] 

or 

(4.3.10)  cos  9  -  -  |  i  lntl  +  j]  . 

Substituting  this  result  in  I  above,  we  obtain  after  some  reduction 

md  +  f) 

(4.3.11)  cos[Xt  +  ln(  +  v)  =  .  _  • 

*  **<1  +  {) 

A  sketch  of  the  expressions  on  either  side  of  this  equation  reveals  the 
locations  of  the  infinity  of  roots.  Some  direct  attempts  at  a  solution 
by  numerical  methods  indicate  that  large  scale  computer  programming  is 
necessary. 

Urns  a  closed  form  approximation  to  the  distribution  function 
Fq(X,T,t)  would  be  helpful.  One  possibility  is  to  fit  a  Gamma  dis¬ 
tribution,  i.e.,  a  density  function  of  the  form 

-ux/  \k-l 

(4.3.12)  r(x>u,k)  -  Ug 

by  fitting  the  first  few  moments  of  this  distribution  to  those  obtainable 
by  differentiating  cpQ(s).  Since  the  first  two  moments  of  r(xjuk)  are 
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(4.3.13) 


k  2  k 
u  «  —  ,  a  =  — * 

0  u  o  7 


and  these  determine  the  particular  Gamma  distribution  completely,  we 
could  set 


(4.3.14) 

and 


[ _ L_.]t  M  £ 

s=o  u 


(4.3.15) 


(CTOT]'32  "  [?TiT]"  =  4  • 

tO  '  S=0  rO  '  8=0  U 


2  2 

If  T  is  large  cpQ(s)  tends  to  X  /(X+s)  ,  that  is,  the  Laplace 
transform  of  a  Gamma  distribution  with  k  =  2  and  X  =  u.  This  is  as 
it  should  be,  for  that  corresponds  to  the  waiting  time  distribution  for 
the  second  input  event  in  the  case  of  a  Poisson  process.  For  large 
enough  values  of  t,  we  are  justified  in  approximating  Fq(X,T,t) 
near  the  origin  by  a  Gamma  distribution  with  k  =  2.  (Some  trial  approx¬ 
imations  easily  show  that  this  parameter  is  not  very  sensitive  to 
variations  in  k.)  Thus  to  a  crude  approximation  we  will  need  only 
u  and  this  parameter  is  found  from 


(4.3.16) 


£=I  + 
u  X 


X(l-e'U) 


or 


(4.3.17) 


2X  l-e  Xt) 
u  =  — '  -Tt 
2-e 


so  that  we  shall  take 
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as  response  function  In  the  case  that  two  photons  within  time  t 
produce  developability.  We  shall  take 


(4,5.19) 


r 


fax(l-e"XT) 


2-e 


-X.T 


J  £ 
1  31 


can 


14'  twice  as  large  a  silver  speck  is  required. 
Formula  (4.3.18)  can  be  integrated  to  yield 


(4.3.20)  Fq(X,T,t)  =  l-e"uT  -uT  e‘uT 

t 

which  indicates  that  the  response  function  of  the  two  photon  photographic 
detector  with  low  intensity  reciprocity  failure  is  identical  in  shape 
(to  this  crude  approximation)  to  that  of  the  two  photon  detector  without 
reciprocity  failure  but  stretched  along  the  X-axis  by  the  factor  X/u. 
Again  to  this  approximation  the  efficiency  of  this  detector  is  related 
to  that  of  the  two  photon  simple  photographic  detector  by  the  relation 

(4.3.21)  <f(X,T,t>  -£(u)  (g)2 

where  <f(u)  is  the  efficiency  of  the  Simple  two  photon  photographic 

detector  evaluated  at  the  abscissa  value  at  which  its  response  curve 

du  2 

has  the  magnitude  Fq(X,T,t) .  The  term  (^)  arises  from  the  fact 

dF  2  dF  2 

that  £(u)  involves  the  expression  (^~)  whereas  we  need  (<jjj^) 
in  the  expression  for  £(X,T,t). 


High  intensity  reciprocity  failure  can  be  handled  in  much  the 
sane  way  as  low  intensity  reciprocity  failure.  In  this  case  we 
postulate  a  dead  time  ,  t,  so  that  only  those  photons  contribute  to 
making  an  R- photon  photographic  detector  respond  which  arrive  at  least 
t  units  apart.  In  this  case  the  probability  that  a  grain  will  become 
developable  in  time  T  is  obtained  by  calculating  the  probability  of 
R  or  more  output  events  from  a  Type  I  counter.  Here  again  we  encounter 
the  apparent  difficulty  that  as  a  first  step  we  need  the  inverse  Laplace 

‘t'Vi 

transform  of  the  R  power  of  the  inter-arrival  time  transform  of 
the  Type  I  counter 

(4.3*22)  [cpo(s)]R  =  e_T8]R 

which  represents  the  waiting  time  distribution  for  the  Rth  output 
event  from  a  Type  I  counter.  Subsequently  we  need  the  integral  from 
0  to  T  which  is  the  probability  that  the  R  event  will  occur 
prior  to  time  T  and  corresponds  to  the  probability  that  a  grain 
exhibiting  high  intensity  reciprocity  failure  will  become  developable. 

However,  in  the  present  case,  we  can  make  use  of  a  device  that 
will  also  be  found  to  be  useful  in  a  problem  to  be  discussed  in  the 
following  chapter.  We  replace  T  by  T'  *  T  -  (R-1)t,  a  contracted 
time  interval.  It  is  easily  shown  that  the  output  process  from  a  Type 
I  counter  is  again  a  Poisson  Process  with  parameter  X.  in  contracted 
time.  Hence  high  intensity  reciprocity  failure  is  equivalent,  under 
the  present  model,  to  shortening  of  the  exposure  time  for  an  R  photon 
simple  photographic  detector  by  an  amount  (R-1)t  --the  accumulated 
dead  time  arising  from  the  first  (R-l)  input  photons  which  are 
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incident  on  a  photographic  grain.  Hence  the  efficiency  of  an  R- photon 
photographic  emulsion  with  high  intensity  reciprocity  failure  and 
time  Is  found  from 

(4.J.23)  H(X,T,t  ,R>  -  H(XT',R) 

and  using 

(4.3.24)  ^  (X,T,Td,R)  - 


we  get 


(4.3.25)  ^(X,T,rd,R) 


X  1  rdfl(X,T,Td,R)i2 

T  ’  H ( X, T, xd,R) [1-R( X,T, td,R) ]  1  dX  1 


which  becomes 


(4.3.26)  <f(^,T,Td,R)  >  .  £(X,T’,R) 

as  can  be  seen  by  multiplying  the  numerator  and  denominator  of 
(4.3.25)  by  T'  =  [T-(R-l)Td]  and  using  (^.3.23)  and  (4. 3. 2k). 


CHAPTHl  V 


COMBINATIONS  OP  COUNTED  ACTING  AS  DETECTORS 

So  far  we  have  been  entirely  concerned  with  detectors  which  are 
employed  to  estimate  the  parameter  X  of  a  Poisson  type  input.  Here 
we  shall  give  some  indication  of  what  can  he  accomplished  in  other 
situations.  First  we  consider  a  series  arrangement  of  a  dead  time  and 
a  simple  coincidence*  counter  which  is  used  as  a  detector  for  discovering 
an  event  consisting  of  two  pulses  with  given  spacing  embedded  in  a 
Poisson  process  of  noise  pulses.  What  is  of  interest  in  such  a  case  is 
the  overall  probability  of  failure,  i.e.,  the  probability  of  not  sensing 
the  twin  signal  pulses  and  the  probability  of  mistaking  two  noise 
pulses  for  the  twin  signal  pulse.  This  can  arise  in  several  physical 
situations,  for  example,  in  a  radar  range  determination  problem,  and 
the  discussion  will  proceed  from  that  point  of  view.  The  second  situa¬ 
tion  in  this  chapter  is  devoted  to  a  detection  situation  in  which  a 
dead  time  and  a  coincidence  counter  are  used  in  parallel.  This  combina¬ 
tion  suggests  itself  in  the  situation  where  the  Poisson  stream  to  be 
observed  consists  of  events  which  have  a  variety  of  energy  levels  and 
it  is  desired  to  estimate  the  intensity  of  that  portion  of  the  stream 
whose  energy  exceeds  a  given  threshold.  This  situation  is  the  discrete 
counter  part  of  a  spectrum  analyzer. 

1.  Dead  Time  and  Coincidence  Counters  in  Series  (Radar  Range 

Determination  Problem) 

A  problem  which  can  be  solved,  at  least  approximately,  by  consider¬ 
ing  a  series  combination  of  dead  time  and  coincidence  counters  acting  on 
a  Poisson  stream,  is  the  following. 
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A  space  vehicle  approaching  the  moon  is  required,  to  fire  its  retro 
rocket  at  some  predetermined  distance  from  the  moon.  To  sense  this 
distance,  pulse  radar  is  used  in  the  following  manner.  The  radar  uses 
two  range  gates.  As  soon  as  a  pulse  is  received  In  the  first  gate  the 
instrument  ceases  to  respond  for  a  time  x^  and  then  fires  the  retro 
rocket  provided  a  signal  pulse  Is  received  in  the  second  range  gate 
during  the  time  interval  x£  following  x^  .  If  no  pulse  is  received 
in  this  interval,  the  first  pulse  is  Judged  to  be  spurious  and  the 
radar  reverts  to  searching  for  the  first  range  indication.  This 
system  can  fail  to  operate  if  (a)  a  spurious  pulse  occurs  Just  prior 
to  the  time  at  which  the  first  turn  indication  would  occur,  thus  im¬ 
mobilizing  and  preventing  it  from  being  in  the  proper  state  to  detect 
the  first  range  signal  and  (b)  two  spurious  pulses  occur  earlier  than 
the  intended  firing  time  spaced  in  such  a  way  as  to  produce  a  pulse  in 
both  range  gates  and  thus  fire  the  retro  rocket  too  soon.  We  are 
interested  in  determining  the  probability  of  both  of  these  types  of 
failures  [p(a)  and  p{b)3  as  a  function  of  the  mean  time  between 
noise  pulses,  l/\,  and  the  total  time  of  operation  of  the  radar 
system,  T. 

The  probability  p(a)  is  very  simply  obtained  and  is  Just  the 

probability  of  a  spurious  event  during  the  time  interval  [T-x^,  T]. 

Since  the  spurious  events  will  be  assumed  to  arrive  in  a  Poisson 

stream  with  parameter  X.,  the  probability  of  no  events  during  this 

-XTj 

time  interval  is  1-e  .  The  probability  p(b)  is  mare  difficult 

to  obtain,  but  an  approximate  value  can  be  obtained  by  the  following 
line  of  reasoning. 
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if 

ftt  system  described  l«  equivalent  to  passing  the  Pol® eon  stream 
first  through  aType  I  counter  vith  dead  time  then  through  a 

hypothetical  time  contractor  which  reduces  the  total  elapsed  time  by 
the  amount  of  dead  time  Mt^  experienced  due  to  the  number,  M,  of 
output  counts  from  the  Type  I  counter  which  have  occurred  during  the 
-elapsed  time.  Ibis  results  in  a  new  Poisson  stream  on  the  Interval 
T'  -  T-»Hpt4  with  parameter 

(5.1.1)  V  - -  . 

E{T' } 


This  stream  is  fed  into  a  coincidence  counter  with  coincidence  time 
constant  and  the  probability  of  one  or  more  coincidences  during 
a  time  span  T' ,  given  that  M  output  pulses  are  obtained  from  the 
first  counter,  can  then  be  obtained  by  integrating  the  inverse  Laplace 
transform  of  the  inter-arrival  time  distribution.  We  then  multiply 
this  probability  by  the  probability  of  getting  M  output  pulses  from 
the  Type  I  counter  and  sum  over  M  to  obtain  the  probability  p(b). 

This  procedure  is  far  from  simple,  however,  and  an  approximation 
to  p(b)  can  be  obtained  by  using  the  expected  number  of  counts  from 
the  Type  I  counter 


(5.1.2) 


SMTH  - 

a 


and  computing  the  probability  of  a  coincidence  for  a  Rjisson  process 
with  parameter  V  over  the  expected  time-contracted  interval 
E(T' )  »  T-E(^}t^.  It  is  interesting  to  note  that 


(5.1.3) 


EtMj) 

E(T') 
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becomes 


(5-1.*) 


XT 

XrTi 

u 

W 


T  - 


ir+T 

CL 


that  is,  the  output  stream  from  the  Type  I  counter  in  the  contracted 
time  interval  T1  is  a  stream  with  the  same  parameter  as  the  input 
process . 

Furthermore,  this  stream  is  a  Poisson  stream  since  the  waiting 
time  distribution  from  the  termination  of  dead-time  (which  corresponds, 
in  contracted  time,  to  the  arrival  of  an  input  event)  to  the  next  input 
event  is  a  negative  exponential  distribution. 

Now  the  expected  number  of  counts  E{M' )  in  a  coincidence  counter 
over  a  time  span  T'  is 

-Xt 

(5-1.5)  E{M'(T'))  =  -C- 

2-e  c 


which  for  Xt  «  1  becomes 


XT'[1-(1-Xt  )]  p 

(5.1.6)  E(M(T'))  =  - -  *•  X  t  T'  . 

2-  (1-Xtc)  c 


Let  P^,, (r)  denote  the  probability  of  exactly  r  coincidences 
in  the  time  T' .  Since  the  expected  number  of  coincidences  is 


oo 

(5.1.7)  E{M(T'))  =  I  rP  (r) 

r=0  1 

and  in  the  present  situation  E{M(T')J  «  1  we  may  make  the  further 
approximation  PT,  (r)  *  0  for  r  >  2  so  that 
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(5.1.8) 


*{M(T'))  -  P  (1)  -  l  P_,(r)  . 

1  r=>l 

The  expression  on  the  right  is  p(b),  the  probability  of  one  or  more 
of  the  kind  of  coincidences  which  cause  premature  firing.  For  this 
approximation,  therefore,  we  have 

(5.1.9)  p(b)  -  X2tcT'  . 

»2 

As  an  example,  let  X  ■  10  counts  per  second  and  T  *  100  seconds, 

t ,  -  t  ■  0.5  seconds  then 
d  c 

-2  2 

(5.1.10)  E{M)  -  «  1 

5xl0°+l 

hence  E(T' )  =  T  (note  that  without  the  coincidence  counter 

we  would  expect  a  noise  pulse  and  hence  premature  firing. )  Since 

E(M(T'))  =  X2t  T  =  5x10“^  ,  the  probability  of  an  early  firing  is 

C  -XTd  -5 

only  about  l/2  percent,  and  since  p(a)  =  1-e  =  5*10  the  over¬ 

all  probability  of  malfunction  due  to  the  incident  BolsBon  noise  is 
p(a)  +  p(b)  *  .01. 

2.  Coincidence  and  Dead  Time  Counters  in  Parallel  -  Energy  Detection 

It  was  remarked  earlier  that  radiation  can  be  considered  as  made 
up  of  discrete  photons  and  that  differences  in  the  frequency  or  wave¬ 
length  of  the  radiation  can  be  associated  with  differences  in  the 
energy  which  the  photons  may  have.  In  some  applications  a  detector  is 
required  to  measure  not  only  the  total  number  of  photons  in  a  stream 
but  also  the  fraction  having  a  given  energy  level.  The  continuous 
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analog  of  this  type  of  a  detector  would  be  a  spectrum  analyzer  which 
measures  the  relative  power  in  a  narrow  band  of  frequencies. 

Here  we  shall  confine  ourselves  to  a  Poisson  stream  in  which  the 
photons  can  have  two  energy  levels  and  Eg.  We  shall  further 
assume  that  the  differences  between  these  types  of  photons  manifests 
itself  in  their  ability  to  produce  electrons  in  some  counter  device; 
the  lower  energy  photons  (E^)  being  able  to  dislodge  only  one  elec¬ 
tron  while  the  higher  energy  photons  are  capable  of  dislodging  at  least 
two  electrons.  The  counters,  whether  dead  time  or  coincidence  time, 
are  conceived  of  as  operating  on  the  electron  stream  produced  by  the 
photons.  We  are  interested  in  observing  both  the  intensity  of  the 
stream  and  the  fractional  intensity  of  the  high  energy  radiations. 

Assuming  that  all  the  electrons  produced  by  a  photon  are  generated 
within  a  time  t  after  the  occurrence  of  the  photon  we  can  proceed  as 
follows: 

A  type  I  dead  time  counter  with  time  constant  t  is  used  to  observe 
the  total  intensity,  X  =  X^  +  Xg,  of  the  overall  Poisson  stream.  With 
this  choice  of  time  constant,  the  multiple  events  produced  by  the  highly 
energetic  photons  will  be  reduced  to  simple  counts  and  the  formulas  of 
Chapter  II  section  2  can  be  used  to  determine  the  efficiency  as  far  as 
the  estimation  of  X  is  concerned. 

In  addition,  a  simple  coincidence  counter  is  introduced  into  the 
stream,  with  coincidence  time  constant  t.  This  counter  will  count  the 
high  energy  photons  (one  count  per  high  energy  photon  since  the  electrons 
are  all  produced  within  a  short  time  t  of  the  occurrence  of  the  photon) 
and  will  also  register  a  count  when  two  low  energy  photons  occur  within 
t  time  units. 
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Let  n^(T)  and  ng(T)  denote  the  number  of  low  and  high  energy 

photons  occurring  during  an  observation  time  T.  An  ideal  device  would 

in  this  case,  produce  output  counts  m*  which,  for  the  overall  X-detector 

would  be  ra|  -  n^+  n2  and  for  the  high  energy,  or  Xg-deteetor,  m*  »  ng. 

'Htie  non- ideal  counts  described  above  would  yield  outputs  m.  and  m 

a  c 

related  to  those  by 

(5-2.1)  md  <  m*  ■  n1  +  n£ 

and 


(5-2.2) 


m  >  m*  =  n, 
c  —  c  I 


Letting  V*(T)  denote  the  ideal  variance  of  the  estimate  of  Xg/X.^  Xg 
and  V(T)  denote  the  variance  of  the  estimate  obtained  from  the  combina¬ 
tion  of  dead  time  and  coincidence  time  counters,  we  can,  as  in  the  case 
of  the  photographic  detector,  define  a  measure  of  detector  efficiency 
by 


(5-2.3) 


=  V*(T)/V(T)  . 


The  variance  of  the  ideal  estimate  of  .■  —  t—  can  be  obtained 

V  2 

from  the  distributions  of  the  ideal  estimates  of  Xg  and  X^  +  Xg. 
We  have 


(5-2.4) 


(^T)2 

ri 


and 

(5-2.5) 


P(\+ 


-(X1+Xg)T  (X1+X2)r+e 

e  (r+s) I 
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or  since  the  two  kinds  of  Poisson  events  are  assumed  to  he  independent 


(5.2,6)  p(X. 


r 

T 


V 


r+Bs 

--5") 


-v 


(VT)1 


-v 


(^T)8 


The  distribution  of  the  ratio  5 


is  consequently 


p(£  =  ~)  =  2  p(Xg=  kr)  p(X  -  ks) 
k*l 

(5-2.7) 

«  -(XT)  (X  T)kr  -(XT)  (X  T)kS 

■  -TET r  e  TESTr 


From  this  formula  we  obtain  the  variance  of  the  ideal  estimate  V*(T) 
by  means  of 


2  2  -  2  2  -(X?T)  <***>*  -(X  T)  (X  T)k8 

V*(T)  "  ril  .1  ^  X*  •  TETi- 


(5.2.8) 


r  oo  oo 


-J 


_  -  -X_T  (X  T)rk  -(XT)  (XT) 

r,r»rr'c<c  IX 

Jl  s?i  r+B  k?i  *  (kr)l  *  (ks)1 


ks 


We  have  not  yet  specified  how  the  actual  estimate  of  Xg/X^+Xg 
is  to  be  obtained  much  less  what  its  variance  will  be. 

To  obtain  an  estimate  of  Xg  once  we  have  obtained  an  estimate 
of  X  =  X^  +  Xg  we  can  proceed  as  follows: 

We  know  that  the  expected  number  of  cotints  from  the  coincidence 
counter  will  be 

(5.2.9)  E{mo)  -  E(na)  +  f(X-Xg,T) 
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where,  using 


(5.2.10) 


v  (X-X  )T[l-e 
f(X-X  ,T)  -  -  2 


■(X-X0)t 


] 


2-e 


-(X-Xg)T 


Using  the  observed  number  of  counts  m  in  place  of  E{m  },  X_T  in 

C  C  fc 

A 

place  of  E(n2)  and  X  in  place  of  X  we  obtain  an  expression 


(5.2.11) 


mc  =  X2T  + 


A  A  -(X-Xp)t 

(X-XgJ'Td-e 


2-e 


-ft-X2h 


A 

which  can  be  solved  for  the  estimate  X^  .  Hence  also  we  obtain  an 

A  .A  A 

estimate  for  the  relative,  intensity  X^/X  ,  since  X  was  obtained 
by  (2.2.7). 

In  principle  we  could  then  obtain  the  variance  of  this  estimate  and 
the  efficiency  of  this  energy  detector.  The  character  of  the  expressions 
makes  it  clear  that  this  involves  extensive  numerical  work  and  perhaps 
Monte  Carlo  simulation. 

In  principle  such  a  combination  of  counters  can  be  used  as  the  dis¬ 
crete  analog  of  a  spectrum  analyzer.  Where  more  than  two  energy  levels 
are  involved  the  coincidence  counter  would  need  to  be  "tuned"  so  that 
only  those  photons  with  an  energy  in  excess  of  a  given  threshold  will 
produce  more  than  one  electron.  By  varying  this  threshold,  E,  the 
proportion  of  photons  having  energy  in  excess  of  E  could  then  be 
observed . 

It  should  be  noted  that  photographic  detectors  have,  to  some  extent, 
the  characteristics  described  above.  High  energy  photons,  for  example 
X-rays,  produce  one  or  more  developable  grains,  since  they  release  a 
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number  of  electron*  in  passing  through  a  photographic  emulsion.  Low 
energy  photons  contribute  to  developability  hut  are  required  in  greater 
numbers  and  over  a  limited  time  span  in  order  to  produce  a  developable 
grain. 
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